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abstract 

This  Semiannual  Technical  Summary  describes  the  Lincoln  Laboratory 
Vela  Uniform  program  for  the  period  1  April  to  30  September  1980. 
Section  I  describes  progress  in  the  development  of  a  prototype  Seismic 
Data  Center.  During  this  report  period,  advances  have  been  made  in 
three  important  subsystems:  the  Seismic  Analysis  Station,  the  Local 
Computer  Network,  and  the  Data  Base  Management  System.  Section  II 
describes  a  series  of  studies  into  the  seismic  processing  algorithms  to 
be  used  at  the  Seismic  Data  Center.  The  results  of  several  investiga¬ 
tions  in  General  Seismology  are  included  in  Sec.  III. 
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SUMMARY 


This  is  the  thirty-third  Semiannual  Technical  Summary  (SATS)  report  describing  the 
activities  of  Lincoln  Laboratory  funded  under  Project  Vela  Uniform.  This  report  covers  the 
period  1  April  to  30  September  1980.  Project  Vela  Uniform  is  a  program  of  research  into  the 
discrimination  between  earthquakes  and  nuclear  explosions  by  seismic  means.  An  important 
recent  emphasis  of  the  project  is  in  the  development  of  data-handling  and  analysis  techniques 
that  may  be  appropriate  for  the  monitoring  of  a  potential  Comprehensive  Test  Ban  Treaty,  pres¬ 
ently  under  negotiation.  During  FY  1979,  Lincoln  Laboratory  completed  a  design  for  a  Seismic 
Data  Center  (SDC)  that  would  satisfy  this  requirement,  and  the  FY  1980  program  consists  of 
beginning  the  development  of  a  prototype  SDC  which  will  test  and  refine  the  various  design  con¬ 
cepts.  Along  with  this  development,  seismic  research  is  continuing,  with  emphasis  on  those 
areas  directly  related  to  the  operations  of  the  SDC. 

To  date,  the  Lincoln  Laboratory  program  has  focused  on  three  of  the  most  important  sub¬ 
systems  of  the  prototype  SDC:  the  Seismic  Analysis  Station  (SAS),  the  Local  Computer  Network, 
and  the  Data  Base  Management  System.  Hardware  acquisition  for  the  SAS  has  been  essentially 
completed,  and  efforts  in  this  area  have  been  focused  on  the  development  of  a  responsive  wave¬ 
form  display  which  will  include  a  smooth  scrolling  capability.  Hardware  for  the  prototype  Data 
Base  Management  System  is  on  order,  but  has  not  yet  been  received.  Efforts  in  this  area  have 
concentrated  on  the  development  of  a  modification  to  the  kernel  of  the  UNIX  operating  system 
that  will  permit  rapid  data  transfer.  The  modified  system,  TUNE,  is  nearing  completion  and 
may  be  utilized  in  the  SAS,  The  Local  Computer  Network  is  being  developed  to  our  specifica¬ 
tions  by  a  subcontractor.  Interactions  with  this  subcontractor  lead  us  to  believe  that  the  design 
of  this  subsystem  is  progressing  well,  and  we  anticipate  delivery  of  the  subsystem  on  schedule. 

Research  into  the  algorithms  that  will  form  the  basis  for  seismic  processing  at  the  SDC  has 
continued  in  several  areas.  On  the  basis  of  their  spectral  characteristics,  it  appears  pi  ^^ible 
to  identify  local  signals  automatically  with  a  high  degree  of  success.  A  set  of  synthetic  arrivals 
has  been  used  to  evaluate  some  existing  automatic  association  algorithms.  There  is  a  substan¬ 
tial  gap  between  the  performance  of  these  algorithms  and  the  much  better  performance  of  a 
human  analyst,  but  future  research  is  expected  to  reduce  this  difference.  Evidence  is  given  that 
it  may  be  necessary  to  routinely  use  master-event  techniques  for  earthquake  location  in  regions 
of  strong  lateral  inhomogeneity,  such  as  island  areas.  We  have  examined  some  of  the  difficul¬ 
ties  involved  in  the  identification  of  short-period  depth  phases,  and  show  that  useful  information 
can  be  obtained  from  long-period  body-wave  records  when  they  are  available. 

Research  in  general  seismology  includes  a  study  of  the  anomalous  Rayleigh  waves  emitted 
by  certain  events  at  the  Eastern  Kazakh  test  site.  The  characteristics  of  these  surface  waves 
can  be  explained  on  the  basis  of  a  combined  explosion/thrust  faulting  mechanism.  A  series  of 
investigations  have  been  carried  out  into  the  direct  inversion  of  waveform  data  to  obtain  infor¬ 
mation  about  the  velocity  and  attenuation  structure  of  the  Earth.  Initial  experiments  have  used 
mantle  waves  of  long  period,  and  results  using  both  synthetic  and  real  data  show  distinct  prom¬ 
ise.  Involved  in  this  study  (and  many  others)  is  the  computation  of  an  average  path-dispersion 
characteristic  for  a  surface-wave  path  that  traverses  several  different  provinces.  Analysis  of 
this  problem  shows  that  computation  of  this  path  average,  at  least  in  certain  cases,  should  not 
be  based  on  simple  geometrical  path  lengths  in  each  region,  but  should  use  some  additional 


parameters  which  may  be  computed. 


M.  A.  Chinnery 
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I.  SEISMIC  DATA  CENTER  DEVELOPMENT 

A.  PROGRESS  IN  SEISMIC  DATA  CENTER  DEVELOPMENT 

During  the  period  covered  by  this  report,  development  has  continued  on  the  various 
subsystems  (described  in  an  earlier  report1)  of  the  prototype  Seismic  Data  Center  (SDC).  The 
main  emphasis  of  the  program  has  been  on  hardware  acquisition,  and  the  development  of  basic 
systems  software. 

Hardware  acquisition  for  the  Seismic  Analysis  Station  (SAS)  has  been  essentially  completed, 
and  efforts  in  this  area  have  been  focused  on  the  development  of  a  responsive  waveform  display 
which  will  include  a  smooth  scrolling  capability. 

Hardware  for  the  prototype  Data  Base  Management  System  is  on  order,  but  has  not  yet  been 
received.  Efforts  in  this  area  have  concentrated  on  the  development  of  a  modification  to  the 
kernel  of  the  UNIX  operating  system  that  will  permit  rapid  data  transfer.  The  modified  system, 
TUNE,  is  nearing  completion,  and  may  be  utilized  in  the  SAS. 

The  Local  Computer  Network  is  being  developed  to  our  specifications  by  a  subcontractor. 
Interactions  with  this  subcontractor  lead  us  to  believe  that  the  design  of  this  subsystem  is  pro¬ 
gressing  well,  and  we  anticipate  delivery  of  the  subsystem  on  schedule. 

Activity  concerning  the  communications  interface  has  been  limited  to  continuing  discussions 
with  DOE  concerning  data  formats  and  data-transmission  methods.  Until  these  questions  are 
resolved,  work  on  this  subsystem  cannot  begin. 

We  are  conscious  of  the  necessity  to  include  adequate  services  and  support  for  the  research 
community.  Efforts  are  currently  under  way  to  formulate  the  requirements  of  this  community. 
Discussions  have  been  held  with  various  individual  seismologists  and  with  the  Panel  on  Data 
Problems  in  Seismology  of  the  National  Academy  of  Sciences. 

M.  A.  Chinnery 
A.  G.  Gann 


B.  SEISMIC  ANALYSIS  STATION  ARCHITECTURE 

The  Seismic  Analysis  Station  (SAS)  is  a  subsystem  of  the  Seismic  Analysis  Center.  It  will 
be  the  focus  of  seismic  interactive  analysis  for  the  system,  providing  facilities  for  the  analyst 
to  scan,  classify,  and  make  measurements  on  large  amounts  of  seismic  waveform  data.  This 
section  will  describe  our  current  thoughts  about  the  SAS  architecture.  Later  sections  summa¬ 
rize  our  work  on  two  prototype  SAS  systems:  Mod  1,  a  system  for  analyzing  parametric  data 
which  is  basically  operational;  and  Mod  2,  a  waveform  display  system  which  is  now  partially 
implemented. 

1.  Basic  Philosophy  of  the  SAS 

Our  goal  is  to  design  a  system  which  is  flexible  enough  so  that  operators,  within  limits,  can 
tailor  the  system  to  their  needs.  We  also  need  a  system  which,  as  we  learn  about  better  methods 
of  seismic  analysis,  we  can  easily  change. 

We  will  try  to  make  the  processes  which  implement  the  functionality  required  very  fine 
grained,  with  simple  interfaces.  There  will  then  be  much  interprocess  communication,  which 
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I  NIX  does  not  support  efficiently.  This  will  require  a  new  service  -  efficient  interprocess 
communication.  We  plan  to  use  the  TUNE  kernel,  now  being  written  for  the  waveform  database, 
as  a  basis  for  Implementation  of  this  function. 

In  addition,  certain  functions,  whose  operations  are  very  obvious  to  the  operators,  must  be 
very  very  fast.  Waveform  scrolling  must  be  fast  and  smooth.  Changes  to  the  parametric  data¬ 
base  must  quickly  appear  on  the  relevant  display.  This  criterion  may,  at  times,  be  in  opposition 
to  the  idea  of  fine-grained  processes.  Any  functions  which  are  to  be  fast  must  be  carefully  de¬ 
signed;  if  they  are  monolithic,  they  still  must  be  logically  internally  fragmented. 

2.  Major  Software  Subsystems 

The  divisions  of  software  components,  as  we  now  see  them,  are: 

User  Interface 

Keyboard  control 
Process  direction 

Interprocess  Communication  and  Control 

Seismic  and  Other  Applications 

Local  Database 

Parametric 

Waveform 

Parametric  Display 

Graphics  Display 
Waveform 
General 

Graphics  supervisor 

"Analog"  Control 
Joystick 
Data  tablet 
Knob  box 

Voice  data  entry  system 

Network  Interface 

UNIX  Modifications 

Interprocess  communication 

Scheduler 

Disk  I/O 

Interface  to  TUNE  kernel 

The  major  thrust  of  this  exercise  was  to  separate  the  functional  components  of  the  system. 
Note  that  if  implementation  is  carried  out  following  the  lines  of  these  divisions,  the  operator's 
input  and  display  devi-es  would  be  logically  separated  from  each  other  and  from  the  computa¬ 
tional  processes  which  use  them.  This  is  similar  to  the  UNIX  philosophy  of  allowing  the  data  to 
flow  between  processes  by  way  of  the  standard  input  and  output.  One  advantage  of  this  scheme 
would  be  that  display  processes  could  be  divorced  from  any  particular  input  control.  For  ex¬ 
ample,  to  indicate  a  particular  line  on  the  parametric  display,  some  operators  might  want  to 
use  the  keyboard  cursors,  while  others,  the  joystick. 

This  also  allows  us  to  experiment  with  the  user  interface.  We  can  easily  provide  different 
faces  to  the  users  with  the  same  internal  functionality. 

There  might  be  exceptions  to  this  philosophy  of  loose  coupling  of  the  processes  and  devices; 
e.g.,  the  waveform  display,  because  of  speed  considerations,  might  have  to  be  more  tightly 
coupled  to  its  input  and  generating  process. 


The  next  two  sections  deal  with  two  key  components  of  the  system.  Interprocess  commu¬ 
nication  and  control  will  shape  the  implementation  of  the  SAS.  1'he  user  interface,  which  is 
intimately  connected  to  process  direction,  shapes  how  the  operator  views  the  SAS. 

3.  Interprocess  Communication  and  Control 

A  necessary  item  in  the  implementation  of  finely  grained  processes  is  the  existence  of  ef¬ 
ficient  interprocess  communication.  We  see  the  need  for  two  types  of  interprocess  communica¬ 
tion:  imperatives,  which  transfer  commands  and  state  information  between  the  processes;  and 
sets  of  data  which  are  shared  among  the  processes  (see  Fig.  1-1). 

4.  Imperatives 

Imperatives  are  the  n.eans  by  which  a  process  directs  another  process  to  perform  a  specific 
task.  They  are  also  the  means  by  which  information  is  carried  from  one  particular  process  to 
another.  Imperatives  are  messages  from  one  process  directed  to  a  specific  other  process;  or 
perhaps,  from  one  functional  unit  to  another,  where  the  specific  process  is  defined  at  another 
level  -  the  user  interface,  for  example. 

These  messages  do  not  invoke  processes  on  their  own.  But,  one  could  use  a  mechanism  by 
which  a  message  to  the  process  control  module  could  start  up  a  process  and  then  forward  a  por¬ 
tion  of  the  message  to  the  new  process. 

Our  thinking  now  is  that  these  messages  are  relatively  short.  When  large  amounts  of  data 
are  to  be  passed,  the  data  sets  will  be  used,  with  the  messages  giving  information  about  them. 

in  addition  to  the  normal  read  and  write  routines  needed  to  implement  this  message  capa¬ 
bility,  we  also  need  a  check  mechanism  (which  may  be  combined  with  the  read)  to  see  if  a  mes¬ 
sage  is  waiting  to  be  read.  This  is  different  from  the  usual  UNIX  method  of  (when  reading) 
waiting  until  there  is  a  message  to  be  read,  reading  it,  and  then  returning  control  back  to  the 
procedure  which  instituted  the  read. 

5.  Named  Data  Sets 

Named  data  sets  are  sets  of  data  which  are  shared  among  the  various  processes.  They  are 
global  "common"  areas  which  are  administered  by  a  process  dedicated  to  their  administration 
and  well  being.  Since  they  are  similar  to  common  areas,  they  must  be  carefully  used  to  avoid 
problems. 

Any  process  may  create,  or  destroy,  a  named  data  set.  All  processes  which  know  the  name, 
may  access  and  change  the  data  within  the  named  data  set.  The  names  are  assigned  at  creation 
time,  and  passed  around  either  by  means  of  imperatives,  or  as  data  within  other  named  data  sets. 

There  will  be  three  types  of  named  data  sets,  each  with  a  different  protocol.  Named  streams 
are  FIFO  data  paths,  similar  to  UNIX  pipes.  They  would  be  faster  than  pipes,  and  allow  a  pro¬ 
cess  to  have  a  multiplicity  of  open  pipes.  A  routine  may  be  needed  to  check  if  the  stream  has  any 
data  in  it. 

Named  parametric  data  sets  are  similar  to  the  Mod  1  data-access  routines.  They  would  be 
record  oriented  and  carry  information  about  their  record  structure.  They  would  have  the  same 
advantages  of  upward  compatibility  as  do  the  Mod  1  data-access  routines. 

Named  waveforms  would  be  the  waveform  transfer  and  storage  mechanism.  A  named  wave¬ 
form  would  contain  a  single  waveform  segment,  along  with  a  header  describing  the  waveform 
structure  -  the  data  structure,  e.g.,  number  of  points,  not  the  seismic  structure. 
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1'hese  named  data  sets  would  be  kept  in  memory.  One  method  of  implementing  them  is  to 
force  each  of  them  to  have  contiguous  memory  locations,  and  to  reserve  one  mapping  register 
for  access  to  the  current  named  data  set.  We  could  have  a  low-priority  process  wtm  h  would  do 
mentor',  ousekeeping  when  the  operator  pauses.  A  reaper  would  have  to  be  implemented,  as 
well  as  a  method  of  backing  least  recently  used  data  sets  off  onto  disk  and  then  back  into  memory 
when  again  required.  This  would  require  some  knowledge  of  which  modules  are  using  which  data 
sets. 

6.  The  Monitor 

A  central  feature  of  the  SAS  organization  is  the  Monitor  module.  The  Monitor  has  two  major 
responsibilities: 

It  dynamically  configures  the  software  of  the  system  at  the  operator's 
and  designers'  commands.  The  Monitor  is  the  method  of  implementing 
the  idea  that  multiple  processes  for  performing  the  same  task  may  exist, 
and  that  the  actual  decision  of  which  process  to  invoke  should  be  left  until 
the  last  moment.  The  operator  can  control  the  assignment  of  processes 
to  functionality,  with  the  full  range  of  features  provided  by  the  UNIX  shell 
assisting  him  in  making  functional  groups  from  fine-grained  processes. 

It  also  serves  as  a  teletype  handler,  with  all  keyboard  data  flowing  through 
it.  Commands  are  interpreted  by  the  Monitor  and  then  passed  on  to  the 
appropriate  process. 

"The  Monitor's  communication  with  the  other  SAS  processes  will  be  implemented  with  the 
multiplexed  I/O  capability  provided  by  Version  7  UNIX.  Basically,  a  multiplexed  I/O  file  is  a 
means  by  which  several  data  streams  can  be  merged  to  appear  as  one.  From  the  multiplexed 
side,  each  message  sent  or  received  from  the  file  will  have  a  channel  identification  associated 
with  it;  from  the  other  side,  eaclt  channel  will  look  like  a  normal  file.  Once  such  a  multiplexed 
file  is  established,  messages  may  be  sent  in  either  direction. 

Whenever  the  Monitor  is  started  (at  the  beginning  of  a  user  session),  it  will  set  up  a  multi¬ 
plexed  file  with  three  channels  for  each  process: 

The  control  channel  will  be  used  to  receive  commands  typed  in  on  the 
teletype  and  information  from  the  sub-processes. 

The  standard  input  channel  will  be  used  to  pass  typed-in  commands  to 
the  sub-processes. 

The  standard  output  channel  will  be  used  to  collect  output  messages  from 
the  monitor  and  the  sub-processes  and  pass  them  on  to  the  teletype. 

Whenever  a  line  is  typed  into  the  teletype,  the  Monitor  will  read  it  and  determine  whether 
the  line  is  a  command  to  the  Monitor  or  one  of  the  running  sub-processes  by  comparing  it  with 
lists  it  maintains,  if  the  command  is  not  found  in  the  command  lists,  it  is  assumed  to  be  a  shell 
command  line,  and  the  Monitor  will  invoke  the  shell  and  pass  it  the  line  for  parsing  (first  setting 
up  the  standard  l/O  channels  as  described).  In  this  way,  all  the  features  of  the  shell,  including 
redirection  of  input  and  output,  will  be  available  under  the  SAS  Monitor. 


U  lien  any  SAS  sub-process  starts  up,  it  is  required  to  perform  several  housekeeping  fune- 
Mons  that  will  provide  the  Monitor  with  necessary  information  about  tiie  sub-pro.  ess.  The  coder 
of  the  sub-proc  ess  will  view  this  linkage  as  a  series  of  simple  subroutine  calls. 

Some  work  has  been  clone  on  a  prototype  Monitor  to  demonstrate  the  feasibility  of  L..ese 

ideas,  but  problems  were  encountered  with  multiplexed  I/O  software  which  is  a  recent  addition 

to  1  MX.  Comparing  notes  with  others  at  a  UNIX  user's  group  meeting,  we  found  that  many 

installations  have  had  trouble  with  this  software.  Various  people  were  working  on  solutions. 

Therefore,  we  decided  to  suspend  work  on  the  prototype  Monitor  and  to  concentrate  on  Mod  2. 

We  assume  that,  bv  the  time  we  are  ready  to  turn  back  to  the  Monitor,  we  will  be  able  to  profit 

from  what  other  UNIX  users  have  learned.  ,,  „  ,,,  , 

H.  S.  Blumberg 

L.  J.  Turek 

C.  MOD  t:  A  PHOTOTYPE  SAS  SYSTEM  FOR  ALPHANUMERIC  DATA 

Mod  1  is  our  first  prototype  for  the  SAS:  a  system  that  deals  only  with  parametric  data  in 
the  form  of  arrival  lists  and  event  lists.  It  was  described  in  detail  in  our  last  SATS.*  Once 
Mod  1  was  operational,  we  felt  that  we  should  lose  no  time  in  going  on  to  other  stages  of  the 
system.  Therefore,  we  decided  that  the  only  further  work  done  on  the  Mod  1  system  would  be 
on  components  that  would  be  carried  over  into  later  systems,  or  any  changes  that  were  absolutely 
necessary  for  Mod  1  to  be  usable.  Tinder  these  guidelines,  some  minor  bugs  were  fixed  and  sub¬ 
stantial  improvements  were  made  to  the  parameter  database  access  routines  described  in  the 
last  SATS. 

In  writing  software  for  Mod  1,  it  had  become  clear  that  the  logical  functions  provided  by  the 
access  routines  were  entirely  appropriate  for  the  types  of  programs  being  written,  and  we  de¬ 
cided  that  these  routines  should  be  carried  on  for  use  in  the  final  SAS.  Therefore,  a  significant 
effort  was  made  to  make  them  more  efficient,  with  a  resulting  speedup  of  approximately  30  per¬ 
cent  and  a  reduction  in  size  of  approximately  25  percent  from  the  initial  version. 

Even  with  this  improvement,  some  of  the  Mod  1  programs  are  still  rather  slow,  but  we  feel 
that  this  can  be  attributed  to  the  limitations  of  the  physical  environment  in  which  Mod  1  must 
operate  -  a  multi-user  UNIX  system  running  on  a  relatively  slow  computer.  With  the  eventual 
conversion  to  the  SAS  configuration  -  a  single-user  system  on  a  larger  computer  -  and  with  some 
additional  indexing  features  yet  to  be  added  to  the  access  routines,  we  feel  that  response  times 
will  be  reduced  to  acceptable  levels. 

In  spite  of  its  current  deficiencies.  Mod  1  is  a  fully  operating  system  which  has  been  used  as 
a  working  environment  for  the  testing  and  improvement  of  such  seismic  programs  as  the  auto¬ 
matic  association  and  location  algorithms. 

1 .  Software  Database 

One  area  in  which  we  hoped  to  gain  experience  during  the  Mod  1  project  was  the  management 
of  a  diverse  software  collection  undergoing  progressive  improvement  by  a  number  of  people. 

Since  the  software  was  to  be  developed  in  stages  and  the  specifications  were  to  change  with 
experience,  it  was  necessary  to  release,  at  appropriate  points,  a  consistent  set  of  software  for 
testing.  A  formal  release  procedure  was  required  so  that  programs  would  match  a  stable  set  of 
documentation  after  a  certain  point.  Programmers  would  generally  like  to  make  "just  one  more 
small  fix  or  enhancement"  before  freezing  a  version  of  a  program.  That  kind  of  flux  makes  co¬ 
ordination  between  modules  and  effective  testing  difficult,  if  not  impossible.  The  diagram 
(Fig.  1-2)  gives  a  generalized  view  of  the  structure  of  the  database  which  we  constructed. 
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The  collection  of  software  which  corresponded  to  a  set  of  documentation  prepared  in  advance 
was  termed  a  "version"  and  given  a  number  (e.g..  Mod  1.1,  Mod  1.2,  .  .  . ).  All  software  for  a 
given  version  was  collected  under  one  branch  of  a  directory  tree.  Within  the  version  directory, 
there  were  three  types  of  subdirectories:  one  for  each  subroutine  library,  one  for  each  pro¬ 
gram,  and  one  for  definition  files.  Within  each  subdirectory,  along  with  the  sources,  we  would 
create  a  command  file  called  "makefile"  which  could  be  run  to  compile  all  the  sources  within  the 
subdirectory.  A  "makefile"  was  also  present  in  the  version  directory  and  contained  commands 
to  run  all  the  subdirectory  command  files.  Thus,  it  was  possible  to  recompile  all  sources  for 
the  version  by  running  the  "makefile"  in  the  version  directory. 

Our  method  of  doing  a  release  consists  of  collecting  all  source  codes  for  an  agreed  upon 
stage  of  the  development  into  one  branch  of  a  directory  tree  structure,  recompiling  the  subroutine 
packages  and  the  programs,  correcting  any  inconsistencies,  making  the  sources  read-only  except 
to  the  release  manager,  and  linking  the  binaries  to  publicly  accessible  names.  Public  names  for 
libraries,  documentation,  and  executable  binaries  carry  a  ".version_number"  suffix  so  that  any 
version  of  a  particular  library,  program,  or  document  can  be  obtained.  In  addition,  the  cur¬ 
rently  released  version  of  each  item  is  available  under  its  name  with  no  version  suffix.  Follow¬ 
ing  the  release  of  a  version,  all  subsequent  software  changes  are  then  incorporated  in  the  next 
numbered  version,  with  the  exception  of  serious  bugs  that  make  it  impossible  to  test  certain 
features  of  the  release.  In  the  case  of  such  a  bug,  the  person  in  charge  of  the  release  super¬ 
vises  the  change  and  recompiles  all  related  software  to  insure  consistency.  It  is  unrealistic  to 
expect  individual  software  authors  to  check  their  own  software  for  consistency  on  a  fixed  sched¬ 
ule  or  to  release  the  software  they  have  developed  in  time  for  a  release  date.  One  individual 
must  be  responsible  for  a  given  software  release  and  must  insure  that  binaries  for  a  given  version 
match  their  corresponding  sources. 

In  doing  a  release,  we  have  found  it  necessary  to  archive  within  the  directory  1  ranch  for  that 
version,  all  subroutine  libraries  and  relevant  files  containing  definitions  used  by  the  programs. 
Should  these  be  changed  for  a  later  version,  it  would  be  impossible  to  recompile  the  original 
sources  (fer  example,  in  the  event  of  a  destroyed  binary)  without  restoring  old  files.  We  also 
found  it  important  to  reference  such  subroutines  and  files  in  such  a  way  that  no  version  numbers 
were  mentioned  in  the  sources.  This  removes  the  necessity  of  changing  references  in  source 
files  when  copying  them  to  the  directory  for  the  next  version.  In  our  archiving  scheme,  the 
directory  structure  for  a  given  version  was  a  subset  of  the  structure  of  subsequent  versions. 

In  general,  we  have  learned  that  it  is  important:  (a)  to  collect  all  files  for  a  particular 
stage  of  a  project  under  one  directory  and  its  subdirectories,  and  to  make  that  collection  as  com¬ 
plete  as  possible  so  that  outside  changes  will  not  destroy  the  consistency  of  the  software  within 
the  branch,  even  though  this  implies  duplication  of  files  in  successive  versions;  (b)  to  automate 
the  compilation  procedure  so  that  it  can  easily  be  duplicated  exactly  at  a  later  date;  (c)  to  keep 
files  as  close  as  possible  to  the  directory  level  of  other  files  that  depend  on  them;  and  (d)  to  avoid 
multiple  copies  of  files  within  a  version,  although  multiple  names  for  the  same  file  can  be  helpful. 

2.  Location  Routine 

The  Mod  1  location  program  uses  a  modified  version  of  the  location  program  written  by 
2 

Julian.  A  "C"  main  program  written  by  Leslie  Turek  provides  the  interface  with  the  rest  of  the 
Mod  1  package.  The  basic  approach  of  the  location  program  remains  intact,  in  that  a  system  of 
linear  equations  is  developed  and  solved  iteratively  by  a  least-squares  method  for  correction 
vectors.  Modifications  are  as  follows: 
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Previously,  the  program  would  weight  each  arrival's  effect  according  to 
its  standard  deviation  and  residual.  This  part  of  the  weighting  scheme 
gave  us  no  trouble.  Unfortunately,  since  this  program  is  to  be  fed  initial 
locations  that  are  potentially  poor,  weighting  by  residual  (especially  dur¬ 
ing  the  first  few  iterations)  had  disastrous  results.  We  found  that  good 
arrivals  would  be  thrown  out,  often  permanently.  Presently,  weighting 
is  dependent  solely  on  the  standard  deviation  of  the  arrival,  which  is  set 
according  to  what  type  of  phase  it  is;  that  is,  typically  impulsive  arrivals 
are  assumed  to  have  been  picked  more  accurately  than  typically  emergent 
arrivals. 

Previously,  we  could  not  process  events  with  less  than  4  arrivals,  as 
this  would  result  in  less  than  0  degrees  of  freedom.  Presently,  if  a 
3 -arrival  event  is  passed  to  the  routine,  it  will  fix  the  depth  at  its  present 
value  or  33  km,  whichever  is  greater,  and  proceed. 

Previously,  the  program  applied  a  damping  factor  to  the  correction  vec¬ 
tors  independent  of  the  number  of  arrivals.  Presently,  the  damping 
factor  is  dependent  on  the  number  of  arrivals,  and  if  there  are  more  than 
20  arrivals  the  corrections  are  applied  full  strength. 

Previously,  if  the  routine  found  that  the  number  of  degrees  of  freedom  was 
less  than  zero,  the  location  process  was  halted.  Presently,  one  of  two 
things  may  happen.  If  the  number  of  degrees  of  freedom  is  equal  to  -1, 
then  the  depth  is  fixed;  if  it  is  less  than  -1,  then  the  window  within  which 
a  residual  must  fall  in  order  for  that  arrival  to  be  considered  usable  is 
increased. 

Presently,  if  the  program  has  used  all  its  allowed  iterations  and  it  still 
has  not  converged  on  a  final  location,  then  the  depth  is  fixed  and  another 
attempt  is  made  to  get  a  stable  location. 

In  order  to  better  simulate  the  Earth,  we  have  added  ellipticity  correc¬ 
tions  to  our  travel-time  calculations. 

The  routine  has  been  translated  into  RATFOR,  which  will  aid  any  further 
attempts  at  modification.  L.  ^  Turek  p>  T_  Crames 

D.  A.  Bach  M.  A.  Tiberio 

D.  MOD  2:  A  PROTOTYPE  SAS  SYSTEM  FOR  WAVEFORM  DATA 

Mod  2  is  going  to  be  our  second  major  milestone  in  the  effort  to  implement  a  complete  SAS. 
Its  purposes  are  to  provide  experience  with  the  Megatek  display  hardware,  provide  a  test  bed 
for  display  software,  and  provide  a  test  bed  for  an  experimental  local  waveform  database. 
Although  it  will  use  much  of  the  parametric  software  written  for  Mod  1,  it  will  not,  at  least  in 
the  early  stages,  interact  in  a  meaningful  way  with  the  Mod  1  database. 

Work  has  been  progressing  along  many  fronts,  all  of  which  will  culminate  in  a  waveform 
display  program,  wdsp.  Among  the  sub-projects  undertaken  and  detailed  below  were:  conver¬ 
sion  to  our  system  of  graphics  routines  obtained  from  Purdue  University's  graphics  project; 
conversion  of  graphics  routines  used  at  our  installation  for  a  Tektronix  storage  display  to  the 
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Megatek  graphics  system;  and  the  writing  of  experimental  programs  to  explore  various  features 
of  tlie  Megatek  which  we  will  need  for  the  waveform  display  program.  In  addition,  the  basic  de¬ 
signs  of  the  waveform  display  program  and  the  waveform  database  have  been  specified.  In  all, 
progress  is  being  made  in  the  three  areas  of  Mod  2:  understanding  the  capabilities  and  finding 
the  limitations  of  the  graphics  display  hardware,  generating  useful  display  software,  and  design¬ 
ing  a  prototype  local  waveform  database. 

1.  Megatek  Graphics  Terminal 

We  recently  purchased  a  Megatek  7000  display  terminal  to  be  used  as  part  of  the  SAS.  The 
Megatek  is  a  high-speed  vector  refresh  display  terminal,  with  hardware  implemented  rotating, 
scaling,  and  translating  of  graphics  segments.  The  refresh  capability  will  be  an  important  fea¬ 
ture,  since  we  will  be  able  to  move  and  scale  waveforms  without  erasing  the  screen,  which  was 
necessary  on  the  Tektronix  4014  terminal. 

We  have  also  obtained  a  Data  Tablet  and  a  Joystick  to  be  used  in  conjunction  with  the  display 
terminal.  The  Data  Tablet  is  a  panel  approximately  1-ft^,  with  a  "pen"  which  is  used  to  identify 
a  location  on  the  tablet.  Megatek  software  can  read  the  tablet  and  translate  the  position  of  the 
pen  on  the  tablet  into  a  position  on  the  screen.  Hence,  we  can  have  a  user  move  the  pen  and 
simultaneously  track  its  position  with  a  cursor  on  the  display  screen,  or  move  a  waveform  on 
the  screen  corresponding  to  where  he  moves  the  pen. 

There  is  also  a  switch  incorporated  in  the  pen  which  is  activated  by  pressing  the  pen  down 
on  the  tablet.  This  switch  will  be  useful,  since  we  can  track  the  pen  with  a  cursor  until  the  user 
presses  down  activating  the  switch.  W'hen  the  switch  is  closed,  we  can  perform  some  function 
such  as  identifying  a  menu  item  on  the  display  terminal,  or  identifying  a  peak  or  trough  on  a 
waveform. 

The  Joystick  is  a  spring-loaded  stick,  normally  in  the  upright  position,  mounted  in  a  small 
box.  When  the  user  moves  the  stick,  Megatek  software  can  detect  the  angle  of  the  stick  and  the 
distance  of  the  stick  from  the  center  position.  On  the  top  of  the  Joystick  is  a  switch  button  which 
can  be  pressed.  Since  the  Joystick  has  features  similar  to  the  Data  Tablet,  it  can  be  used  in  a 
similar  manner.  In  the  future,  we  plan  to  acquire  a  Function  Switches  and  Control  Dials  module 
(FS CD).  The  FS CD  is  a  general-purpose  device  with  eight  rotatable  knobs  and  sixteen  switches 
with  lights.  This  device,  under  software  control,  will  be  used  to  move  waveforms  within  a  fixed 
window,  to  identify  points  on  a  waveform,  and  for  general-purpose  communication  between  the 
user  and  the  software. 

2.  Purdue  Routines 

When  we  acquired  the  Megatek  7000,  we  also  obtained  a  library  of  graphics  subroutines 
written  for  the  Megatek  from  Purdue  University.  Included  in  the  Purdue  routines  are  functions 
to  plot  lines,  move  cursors  and  segments,  read  the  Data  Tablet  and  Joystick,  and  display  text 
strings.  These  routines  are  written  entirely  in  "C".  The  Purdue  routines  allowed  us  to  imme¬ 
diately  test  the  Megatek,  and  to  develop  several  demonstration  programs.  These  routines  were 
also  used  to  convert  Graphpac,  a  set  of  graphics  routines  developed  for  the  Tektronix  4014,  to 
run  on  the  Megatek.  Besides  supplying  graphics  routines,  Purdue  University  also  sent  us  a 
device  driver  for  the  Metatek.  We  modified  it  to  conform  with  our  hardware  and  version  of 
UNIX,  and  then  installed  the  graphics  routines. 
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3.  Graphpac  Conversion  to  Megatek 

About  one  year  ago,  we  developed  a  set  of  subroutines  and  functions  to  perfor  m  graphics  on 
the  Tektronix  4014  graphics  terminal,  called  Graphpac.  These  routines  were  described  briefly 
in  the  previous  SA'TS.1  Since  they  were  written,  several  seismologists  as  well  as  computer  per¬ 
sonnel  have  used  Graphpac  to  build  displays  containing  axes,  waveforms,  and  text.  With  many 
programs  running  with  Graphpac  on  the  'Tektronix,  we  thought  it  would  be  useful  to  convert 
Graphpac  so  that  it  could  be  used  on  the  Megatek.  An  advantage  of  having  Graphpac  on  the 
Megatek  is  that  Graphpac  routines  can  be  called  from  Fortran  (as  well  as  "C"),  which  is  the 
language  many  seismologists  prefer  to  use.  Because  there  are  only  a  few  routines  in  Graphpac 
which  actually  output  data  to  the  Tektronix,  the  conversion  went  smoothly.  The  Purdue  graphics 
routines  for  the  Megatek  were  used  as  an  interface  between  Graphpac  and  the  Megatek  7000. 
During  the  conversion,  there  were  several  obstacles  which  we  overcame;  these  included  differ¬ 
ent  screen  sizes  between  the  Tektronix  and  the  Megatek,  different  methods  of  addressing  the 
screen,  and  simulating  the  Tektronix  crosshairs  with  the  Megatek  Data  Tablet.  All  software 
which  runs  with  Graphpac  on  the  Tektronix  can  now  be  run  on  the  Megatek  simply  by  recompiling 
the  program  with  the  converted  routines. 

4.  Initial  Experiments  with  the  Megatek  Graphics  Terminal 

A  number  of  programs  were  written  to  test  various  aspects  of  the  Megatek  display  hardware, 
and  its  interreaction  with  the  DEC  PDP-11  and  the  UNIX  operating  system.  Among  the  programs 
written  was  one  to  scroll  sine  waves  at  a  speed  and  direction  determined  by  the  Joystick. 

This  was  a  very  useful  program.  It  revealed  two  problems  at  an  early  stage  which  we  will 
have  to  deal  with.  First,  and  most  easily  disposed  of,  was  that  the  Joystick  does  not  transmit 
a  fine  enough  indication  of  its  position.  This  can  be  solved  in  two  ways:  first,  by  replacing  the 
Joystick  with  a  more  sensitive  one,  or  by  using  a  different  control  method,  e.g.,  a  knob.  We 
expect  to  interface  a  knob  control  to  the  system  at  an  early  stage. 

The  second  problem  is  not  as  easy  to  resolve.  We  found  that  the  UNIX  operating  system, 
even  when  supporting  a  single  user,  was  not  consistent  enough  to  present  a  pleasing  moving  dis¬ 
play.  The  rate  of  movement  was  too  variable  for  an  operator's  comfort.  To  use  it  as  it  stands 
would  result  in  a  decrease  in  operator  productivity. 

The  crux  of  the  problem  is  that  the  display  computer  must  be  updated  in  a  very  consistent 
manner  to  support  waveform  scrolling.  Were  it  not  for  the  scrolling  requirement,  the  existing 
software  would  be  more  than  adequate. 

There  are  a  number  of  methods  by  which  we  believe  we  can  solve  this  problem.  We  can 
write  a  "smart  display  driver"  which,  from  within  the  kernel,  controls  waveform  scrolling.  We 
can  use  the  TUNE  kernel  which  is  being  written  for  the  waveform  database  subsystem.  Or,  we 
can  use  another  operating  system  which  is  more  amenable  to  real-time  programming. 

The  last  option,  a  change  in  the  operating  system,  seems  the  least  attractive.  It  would 
require  a  high  overhead  to  support  an  operating  system  which  would  be  unique  within  the  Applied 
Seismology  Group.  Also,  the  available  operating  systems  are  not  as  supportive  of  program 
development  as  is  the  the  UNIX  operating  system  which  we  are  now  using. 

The  best  option  seems  to  be  to  use  TUNE  with  co-supervisors  of  UNIX  and  s  special  graphics 
mini-supervisor.  Tl'NE  would  handle  the  communication  between  them  as  well  as  providing  the 
real-time  primitives  which  are  required.  This  would  fit  into  our  plans,  as  we  had  expected  to 
use  TUNE  at  a  later  date  for  efficient  interprocess  communication. 
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Since  I't  NK  with  a  (  MX  supervisor  is  not  yet  ready,  we  have  implemented  a  "semi- 
intelligent"  graphics  driver  as  a  temporary  measure.  This  graphics  unit  driver  can  queue  re¬ 
quests  to  change  the  display  list,  and  will  transmit  one  for  each  sync  signal  received  from  the 
graphics  display  unit.  Each  change  request  consists  of  a  block  i  f  changes,  at  different  places 
within  the  display  list.  I’ he  sync  signal  is  transmitted  by  the  graphics  display  unit  after  it  com¬ 
pletes  each  scan  of  the  display  list.  The  driver  then  changes  the  display  list  between  the  time 
the  scan  ends,  and  the  next  one  begins.  Thus  real-time  waveform  scrolling  can  be  implemented 
as  synchronized  changes  to  the  display  list. 

Also  implemented  is  a  command  to  clear  the  queue.  This  is  necessary  to  stop  or  change  the 
speed  of  the  waveform  scrolling  in  a  timely  fashion.  This  graphics  unit  driver  will  probably 
enable  us  to  produce  a  Mod  2  waveform  display  program  which  will  be  responsive  and  reasonably 
pleasant  to  use. 

5.  The  Waveform  Display  Program 

The  culmination  of  Mod  2  will  be  a  working  waveform  display  program.  This  program  will 
display  up  to  twelve  waveforms  on  a  screen  which  has  been  divided  into  a  number  of  horizontal 
bands.  This  program  is  not  meant  to  be  used,  at  least  initially,  for  seismic  analysis.  The  com¬ 
mands  of  the  program  have  been  specified,  the  design  process  is  now  proceeding,  and  some  of 
the  service  routines  have  been  written. 

The  commands  consist  of  one  or  more  waveform  numbers,  a  command  letter,  and  a  variable 
number  of  arguments.  A  command  syntax  analysis  routine  has  been  written  and  tested  which 
reads  the  command  from  the  standard  input  device,  analyzes  it,  and  constructs  a  structure  which 
contains  the  information.  It  also  does  error  checking  and  only  allows  valid  commands. 

Among  the  commands  to  be  implemented  are: 

d  Display  a  named  waveform. 

m  Horizontally  move  a  waveform,  or  group  of  waveforms, 
to  a  specific  time,  or  relative  to  the  present  time. 

o  Overlay  two  waveforms. 

s  Horizontally  scroll  a  waveform  or  a  group  of  waveforms. 

k  Mark  a  waveform;  i.e.,  attach  a  marker  to  the  waveform  data. 

This  marker  contains  an  ASCII  string  to  identify  it. 

r  Rescale. 

e  Exit. 

This  set  of  commands  will  enable  us  to  test  the  concepts  of  the  graphics  programming  we 
proposed.  It  will  also  gi-,e  us  some  early  experience  with  person/machine  interfaces.  The 
primitives  needed  to  create  this  program  will  be  used  later  to  implement  the  final  SAS. 

The  early  design  showed  that  we  need  a  block  allocation  method  for  the  graphics  unit  mem¬ 
ory.  A  group  of  memory  management  routines  to  perform  this  task  have  been  written  and  tested. 
They  perform  in  much  the  same  way  as  the  allocation  scheme  for  the  UNIX  memory  allocation 
of  the  PDP-11. 
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t>.  Graphics  Display  l  nit  Language  Assembler  and  Linker 

Another  requirement  which  was  recognized  early  in  the  design  process  of  the  waveform 
display  program  was  a  need  for  an  easy  method  of  specifying  standard  pictures.  There  is  a 
need  for  a  set  of  static  sketches  static  in  the  sense  that  the  specifications  of  the  picture  do 
not  change  as  opposed  to  the  fact  that  the  picture  remains  in  one  position  or  is  even  displayed 
on  the  screen.  1'hese  static  sketches  are  analogous  to  library  subroutines.  They  can  be  com¬ 
bined  with  other  sub-pictures  to  giv  e  a  complete  display.  Examples  of  these  are  the  graphics 
associated  with  waveform  marks  and  axes. 

We  have  decided  to  build  a  Megatek  Assembler  and  Linker.  The  assembler  would  accept  a 
simple  Megatek  assembly  language,  a  symbolic  representation  of  its  display  list  language.  The 
assembler  would  then  analyze  this  input  and  produce  a  file  which  is  the  octal  representation  of 
the  input  display  list.  1'he  linking  loader  would  be  a  subroutine  callable  from  the  waveform  dis¬ 
play  program.  It  would  read  the  output  file  of  the  assembler,  load  the  display  lists  into  Megatek 
memory,  and  return  to  the  calling  program  a  table  of  addresses  where  the  various  display  lists 
could  be  found. 

1'he  Graphics  Assembler  is  comprised  of  five  modules: 

a  lexical  analyzer, 
a  simple  parser, 
a  preliminary  assembler, 
a  runtime  assembler,  and 
a  linking  loader. 

1'he  lexical  analyzer  tokenizes  the  assembler  source  input  file,  dividing  it  into  meaningful  pieces 
for  the  parser,  [  he  parser  performs  the  actual  interpretation  of  the  tokenized  assembler  source. 
It  issues  instructions  to  the  preliminary  assembler  as  to  how  it  should  produce  a  character  graph¬ 
ics  file. 

The  runtime  assembler  reads  in  these  character  files  and  produces  binary  images  with  pos¬ 
sible  external  references  which  the  linking  loader  must  resolve.  The  loader  also  moves  the  re¬ 
sulting  binary  graphics  "object"  to  the  graphics  device. 

The  assembler  source  format  is  not  dependent  on  fixed  columnation;  spacing  is  insignificant 
except  to  separate  the  tokens.  It  does  not  matter  which  columns  of  the  line  a  particular  field 
occupies.  Each  source  line  comprises  source  for  a  single  instruction  and  consists  of  an  optional 
label,  an  instruction  and  its  arguments,  and  an  optional  comment.  In  general,  a  single  mne¬ 
monic  assembly  language  instruction  will  assemble  to  any  one  of  a  variety  of  graphics  processor 
machine  instructions,  depending  on  what  the  arguments  to  the  instruction  are. 

SAMPLE  INSTRUCTION  SET 

VEC  Draws  a  vector,  arguments  indicate  the  end  point  of  the  vector. 

POINT  Plots  a  single  point  on  the  screen,  arguments  indicate  the 
position  of  the  point. 

JUMP  Proceed  to  specified  set  of  graphics  commands. 

JSUB  Execute  specified  set  of  graphics  commands  and  then  resume 

processing  at  this  point. 
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ORIGIN  Set  the  origin  of  the  display  to  the  point  specified. 

I'KX  I'  Display  the  text  string  argument. 

1’MK  Set  tlie  elements  of  the  hardware  vector  transformation  matrix. 

7.  Local  Waveform  Database  Design 

We  have  developed  a  design  for  a  local  waveform  database  to  reside  on  the  SAS  host  com¬ 
puter  and  serve  as  input  to  the  waveform  display  program  and  other  seismic  processing  pro¬ 
grams.  Hat  h  database  will  be  named,  so  that  more  than  one  waveform  database  can  be  in  exis¬ 
tence  at  a  given  time.  In  the  final  implementation  of  a  data  center,  a  local  waveform  database 
will  be  created  when  the  user  initiates  a  data  request  from  the  SAS  to  the  waveform  database 
subsystem.  While  we  are  still  working  with  an  isolated  SAS,  we  have  programs  to  convert  data 
from  the  -gi/.g  format  currently  in  use  on  the  research  1  NIX  system  and  will  provide  programs 
to  read  in  some  types  of  waveform  data  directly  from  magnetic  tape. 

The  waveform  database  designed  for  Mod  2  will  be  analogous  to  the  waveform  database  cur¬ 
rently  in  use  on  our  research  UNIX  system,  with  some  modifications  to  eliminate  problems 
which  have  been  noted  in  the  past.  Logically,  a  wavefor  m  database  consists  of  a  numbered 
series  of  waveform  segments,  with  associated  descr  iptive  infor  mation,  sample  markers,  and 
dropout  flags. 

The  description  file,  which  in  the  old  format  was  a  fixed-field  alphanumeric  file,  will  now 
be  a  parameter  file  which  will  allow  the  use  of  the  parameter  file  access  routines  developed 
under  Mod  1.  Parameters  comprising  the  waveform  dest  ription  file  are: 


wfnn 

int 

waveform  number 

station 

char(8) 

station  code 

channel 

char(8 ) 

channel  code 

stime 

tim_t 

start  time  of  waveform  segment 

etime 

tim_t 

end  time  of  waveform  segment 

spsamp 

int 

(■sampling  frequency 

spsecs 

int 

(.samp1  samples  per  ’spsecS’  seconds 

cal 

float 

calibration 

format 

char 

code  for  data  format  =  32 -bit  integer 

Uach  waveform  data  segment  will  be  a  binary  file  containing  a  continuous  stream  of  samples 
in  the  designated  format.  Data  dropouts  (data  samples  not  received)  will  appear  with  the  value 
of  0. 

Markers  are  a  means  by  which  an  operator-assigned  name  is  associated  with  a  particular 
data  sample  of  a  waveform.  They  can  indicate  phase  onset  times,  peaks  and  troughs,  or  any¬ 
thing  else  required.  Markers  can  be  displayed  or  can  be  used  as  an  input  to  seismic  processing 
programs.  In  the  new  waveform  database,  there  will  be  a  separate  marker  file  for  each  wave¬ 
form  segment.  It  will  be  in  the  form  of  a  parameter  file,  with  only  two  entries  in  each  record: 
the  name  of  the  marker  and  the  sample  number  it  corresponds  to. 

The  dropout  file  is  a  new  addition.  It  will  be  a  parameter  file,  with  each  record  describing 
a  data  dropout  or  glitch  in  the  associated  waveform.  Parameters  will  be  the  starting  and  ending 
sample  numbers  and  a  code  for  the  type  of  problem  (data  never  mceived,  data  judged  to  be  bad 
by  analyst,  etc.). 
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For  purposes  of  fast  scrolling  by  the  display  program,  it  may  be  necessary  to  have  an  addi¬ 
tional  type  of  file  which  incorporates  both  the  marker  and  dropout  information  in  a  form  that  can 
be  accessed  more  rapidly.  1'his  may  not  become  a  permanent  part  of  the  database,  but  may  be 
created  only  for  use  by  the  display  program.  We  call  this  file  the  "bit-string  file"  because  it 
will  simply  contain  one  or  two  bits  to  indicate,  for  each  sample  in  the  waveform,  whether  there 
is  some  mention  of  it  in  the  corresponding  marker  or  dropout  file.  During  fast  scrolling,  this 
file  can  be  read  along  with  the  new  data  with  minimal  head  movement,  and  points  of  significance 
can  be  roughly  marked  as  they  flash  by.  When  the  scrolling  stops,  the  display  program  can 
then  take  the  time  to  retrieve  the  detailed  information  from  the  marker  and  dropout  files  and  to 
display  it  fully. 

At  our  current  stage,  we  need  very  little  special-purpose  software  to  work  with  these  wave¬ 
form  databases.  Most  of  the  files  can  be  read  and  written  using  either  the  parameter  access 
routines  or  regular  binary  reads  and  writes.  We  plan  to  provide  some  simple  tools,  such  as  a 
reference  file  containing  structure  definitions  for  the  various  file  types,  subroutines  to  convert 
from  sample  times  to  sample  numbers  and  back,  and  some  simple  subroutines  to  create  and 
unpack  the  bit-string  files. 

If  this  waveform  database  format  proves  adequate  for  our  needs  in  Mod  2,  we  will  then  de¬ 
sign  a  higher-level  software  interface  oriented  toward  the  needs  of  the  seismic  researcher. 

Such  an  interface  would  include  the  concept  of  a  waveform  segment  (with  all  its  associated  in¬ 
formation)  as  a  data  object  that  can  be  displayed  and  manipulated. 

R.  S.  Blumberg  M.  K.  Nahabedian 

L.  J.  Turek  P.  T.  frames 

E.  DATABASE  SUBSYSTEM:  DESIGN  OF  THE  TUNE  KERNEL 

Development  efforts  on  the  waveform  database  have  focused  on  the  design  and  implementation 
of  the  TUNE  kernel.  Presented  below  is  a  series  of  highlights  of  the  kernel  design  together  with 
information  regarding  the  current  state  of  the  implementation. 

1.  Goals 

The  major  goal  of  TUNE  is  to  provide  a  network  environment  on  a  single-processor  system. 
Within  this  environment,  sets  of  processes  may  operate  as  individual  isolated  host  units  just  as 
they  would  on  their  own  host.  Different  types  of  applications  should  be  provided  proper  support 
so  that  each  may  be  written  in  such  a  fashion  that  it  may  be  transferred  to  its  own  host  with  mini¬ 
mal  change.  The  emulation  environment  should  provide  the  ability  to  test  systems  under  various 
host  and  network  failure  scenarios,  and  also  provide  a  performance  testing  capability  under  a 
wide  variety  of  configuration  assumptions. 

2.  General  Structure 

The  kernel  provides  to  each  supervisor  a  virtual  machine  similar  to,  but  not  identical  with, 
a  PDP-11.  The  differences  arise  both  from  idiosyncrasies  within  the  PDP-11  architecture  and 
from  efforts  to  provide  a  more  felicitous  interface  to  the  supervisor  with  the  goal  of  eliminating 
overhead  and  supervisor  complexity,  without  eliminating  useful  flexibility.  Reduction  of  over¬ 
head  is  particularly  important,  since  TUNE  is  to  be  used  in  environments  requiring  high  I/O 
throughput  which  cannot  be  achieved  while  doing  virtual  machine  simulation  of  the  PDP-11 
architecture. 
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r(  .\i;  provides  kernel  support  for  multiple-user  processes  under  control  of  a  supervisor. 

The  priorities  of  these  processes  and  tile  address  spaces  associated  with  them  are  controlled  by 
the  supervisor,  this  control  being  exercised  by  means  of  explicit  calls  for  kernel  services  which 
are  made  via  the  DM  1'  instruction.  The  user  tasks  may  request  supervisor  services  by  means  of 
the  1KAT  instruction,  which  results  in  an  interrupt  fielded  by  the  supervisor. 

Besides  the  KM  1’  instruction,  there  are  two  other  channels  of  supervisor-kernel  communica¬ 
tion.  A  shared  memory  area,  the  supervisor  interface  area,  is  mapped  into  a  reserved  page  of 
both  the  supervisor  and  kernel  address  spaces  and  serves  as  a  medium  of  low-overhead  commu¬ 
nication.  Supervisor  interrupts  are  used  by  the  kernel  to  notify  the  supervisor’  of  an  external 
event  by  diverting  the  flow  of  control  within  a  supervisor  to  a  supervisor-specified  routine  which 
then  may  deal  with  the  event. 

1.  !  'O  Support 

I'he  n  Mi  kernel  provides  I/O  facilities  which  aim  to  remove  much  of  the  complexity  of 
interrupt  handling  and  device  control  from  the  domain  of  the  supervisor,  while  retaining  as  much 
flexibility  as  possible  for  the  implementation  of  diverse  supervisors  with  diverse  purposes. 

1'hree  different  types  of  I/O  are  distinguished,  each  handled  in  a  different  fashion.  The  three 
types  are  terminal  I/O,  block  I/O,  and  local  network  I/O. 

The  terminal  handling  code  within  the  kernel  is  responsible  for  the  handling  of  character- 
lev  el  interrupts,  only  interrupting  the  supervisor  when  truly  necessary.  Individual  terminal 
interface  modules  provide  a  number  of  different  styles  of  terminal  handling  to  the  supervisors. 

The  supervisors  communicate  with  the  terminal  interface  by  setting  and  interrogating  terminal 
attributes  defined  by  the  individual  interfaces.  The  supervisors1  communication  with  the  termi¬ 
nal  interfaces  is  mediated  by  an  KMT  service  routine  which  provides  for  the  setting  and  interroga¬ 
tion  of  several  attributes  as  an  indivisible  operation. 

The  block  I/O  interface  provides  the  supervisor  witli  a  simple  means  of  controlling  disks 
and  tapes,  which  isolates  the  device-specific  handling  within  a  kernel  driver  routine  yet  provides 
full  support  for  the  asynchronous  operations  required  by  the  waveform  database.  Each  virtual 
device  owned  by  a  supervisor  is  represented  by  a  device  image  located  within  the  supervisor 
interface  area.  Space  is  available  within  the  interface  area  for  the  supervisors  to  construct  I/O 
request  blocks  which  specify  I/O  operations  to  be  performed.  A  queue  of  such  request  blocks 
represents  the  I/O  operations  to  be  performed  on  the  virtual  device.  The  supervisor  may  in¬ 
voke  an  I/O  operation  by  placing  an  I/O  request  block  at  an  appropriate  point  within  the  queue  and 
notifying  the  kernel  of  its  presence  bv  an  KMT  instruction.  The  supervisor  is  free  to  choose  in 
what  order  the  I/O  operations  will  be  done.  When  each  I/O  operation  is  completed,  the  next  I/O 
request  is  taken  from  the  associated  queue.  The  supervisor  is  notified  of  the  completion  an  I/O 
operation  by  a  supervisor  interrupt  or  by  waking  a  specified  process. 

The  interface  to  the  local  network,  both  simulated  and  real,  has  not  been  fuli.v  specified. 

The  following  considerations  will  govern  the  shape  of  this  interface.  It  should  provide  the  same 
flexibility  and  facilities  for  asvm  hronism  as  the  block  1  A)  interface.  It  should  provide  trans¬ 
parent  operation  irrespective  of  whether  it  is  communicating  with  a  supervisor  on  the  same  or  a 
different  host,  while  providing  verv  high  performance  when  used  between  two  supervisors  on  the 
same  host  in  order  to  limit  the  performance  penalty  which  must  be  absorbed  in  handling  small 
configurations  in  this  fashion. 
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4.  Virtual  Clocks 

Both  the  block  I/O  and  local  network  interfaces  provide  for  the  maintenance  of  supervisor 
virtual  clocks,  which  predict  the  time  that  the  supervisor's  execution  would  require  on  its  own 
host  under  the  assumptions  of  the  simulation.  Virtual  clocks  are  maintained  by  determining  the 
time  that  the  primitive  operations  of  a  given  supervisor  will  take  under  the  given  circumstances 
and  using  these  individual  predictions  to  govern  the  evolution  of  the  simulation.  This  requires 
that  interrupts  and  other  items  of  kernel-supervisor  communication  be  assigned  virtual  times 
for  their  performance.  This  kernel  uses  these  times  to  maintain  a  queue  of  events  to  be  per¬ 
formed  ordered  by  virtual  time,  and  synchronizes  processor  access  for  the  given  supervisor's 
processes  with  these  events. 

In  addition  to  the  intra-supervisor  virtual  time  synchronization  required  by  block  I/O,  the 
simulation  of  local  network  operations  requires  inter-supervisor  synchronization,  to  prevent  the 
virtual  clocks  of  the  several  supervisors  from  diverging  and  thereby  invalidating  the  simulation 
as  a  whole.  TI  NE  provides  for  a  specifiable  window  which  specifies  the  maximum  permissible 
divergence  among  all  virtual  clocks  within  a  given  host. 

5.  Implementation  Status 

Considerable  code  has  been  written  for  the  TUNE  kernel.  Some  of  this  has  been  tested.  The 
tests  to  date  have  been  rather  primitive.  Now  that  the  basic  facilities  of  TUNE  are  working,  the 
test  supervisors  can  assume  greater  sophistication  and  test  the  kernel  more  rigorously.  The 
following  areas  have  been  coded  and  tested  fairly  thoroughly:  the  EMT  interface,  system  initial¬ 
ization,  and  address  space  management.  The  following  areas  have  been  coded  and  have  been 
tested  to  some  degree:  process  management,  virtual  clock  management,  and  terminal  handling. 
Two  areas  remain  to  be  coded:  block  I/O  and  local  network  support.  Work  on  these  will  proceed 
once  testing  of  terminal  handling  and  process  management  is  complete. 

D.  B.  Noveck 


F.  LOCAL  COMPUTER  NETWORK 

The  local  computer  network  is  needed  to  interconnect  the  computers  of  the  Seismic  Data 
Center  (SDC)  to  provide  an  effective  and  efficient  means  of  transmitting  data  and  control  infor¬ 
mation  between  computers.  This  computer  network  must  serve  to  interconnect  all  the  comput¬ 
ers  which  will  make  up  the  SDC.  It  must  operate  very  reliably  and  with  a  high  intercomputer 
data  rate. 

The  proposals  for  development  of  the  local  computer  network  were  received  and  evaluated 
during  this  report  period.  Four  proposals  were  received  from  industry  offering  several  ap¬ 
proaches  to  meeting  the  requirements  of  the  local  computer  network.  The  winning  proposal  was 
submitted  by  Sytek,  Inc.  A  contract  for  the  local  computer  network  was  initiated  2  July  1980, 
and  the  first  design  review  was  held  in  August.  The  projected  completion  of  the  system  is 
30  June  1981. 

'  he  Sytek  plan  is  to  develop  a  network  front-end  processor  based  on  a  microcomputer.  This 
processor  will  implement  the  protocols  needed  to  transfer  data  between  hosts  over  a  2-Mbps 
shared  cable  using  CATV  (cable  access  television)  technology  for  data  transmission.  The  data 
transmission  will  use  the  "mid-split"  scheme  with  a  standard  CATV  channel  pair  -  one  fre¬ 
quency  to  transmit,  and  a  related  frequency  to  receive.  At  the  "head-end"  (CATV  terminology), 
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a  frequency  translator  will  convert  all  received  frequencies  up  about  150  MHz,  and  retransmit 
them  on  the  same  cable,  l'lie  separation  of  inbound  and  outbound  signals  is  by  standard  CATV 
filters,  l'lie  mid-split  frequency  translator  is  also  a  standard  CATV  product.  1’he  use  of  CATV 
technology  allows  the  extension  of  the  network  to  the  order  of  kilometers  of  extent,  and  the  sig¬ 
naling  speed  to  be  increased  to  the  order  of  several  tens  of  megabits  per  second.  The  cable  TV 
industry  has  developed  highly  reliable  amplifiers,  frequency  translators,  cable  tapping  hardware, 
etc.,  which  are  now  readily  available  for  use  in  data  systems.  The  use  of  CATV  hardware  with 
\  III'  (very  high  frequency)  modems  for  local  computer  networks  was  pioneered  by  the  MI  TRE 
Corporation.^ 

l'lie  Sytek  approach  includes  the  use  of  a  network  front-end  processor  implemented  by 
Intel  8086  and  8080  chips  to  perform  most  of  the  logic  needed  to  communicate  over  the  network. 
The  808t>  will  implement  the  control  of  the  protocols  and  the  8089  will  control  the  DMA  move¬ 
ment  of  the  data  over  the  Cnibus  to  the  front-end  memory  and  from  the  front-end  memory  to 
and  from  the  network  transceiver. 

The  protocols  to  be  implemented  include  a  link  level  protocol  which  is  a  listen-while-talk 
protocol  similar  to  that  implemented  by  MI  TRE  in  their  MITRENET  system.  The  key  element 
in  this  protocol  is  the  transmission  of  a  preamble,  and  the  detection  of  the  correct  reception  of 
the  preamble  for  sufficient  time  to  allow  all  nodes  to  detect  that  the  system  is  in  use.  If  a  col¬ 
lision  is  detected,  all  transmitting  hosts  stop  and  wait  for  a  small  variable  period  before  trying 
again. 

The  next  protocol  in  the  hierarchy  uses  the  link  protocol  and  provides  the  basic  transmission 
of  data  from  one  node  to  another.  This  is  called  a  transport  level  protocol  and  provides  unre¬ 
liable  (i.e.,  no  acknowledgments)  communication  between  network  front-end  processors.  This 
protocol  is  used  by  the  network  front-end  processors  to  implement  a  reliable  datagram  protocol, 
called  the  transaction  protocol,  and  a  reliable  stream  protocol  used  to  transmit  large  blocks  or 
streams  of  data  from  one  host  to  another. 

The  host  interface  is  implemented  as  a  set  of  Unibus  registers  through  which  command  mes¬ 
sages  are  passed.  The  command  messages  are  composed  in  the  host  memory  and  the  address 
of  the  command  is  passed  to  the  front-end  processor.  The  front  end  then  schedules  the  operation 
and  interrupts  the  host  when  the  operation  is  successfully  completed  or  has  failed.  The  host  then 
can  take  corrective  action  for  failures,  and  schedule  further  operations  for  successes. 

The  development  of  the  local  network  is  proceeding  on  schedule.  The  delivery  of  hardware 
is  scheduled  for  next  spring,  with  a  system  test  period  scheduled  to  end  in  June  1981  with  the 
completion  of  the  contract.  .  „  _ 


G.  SYSTEM  DOCUMENTATION:  ADAPTATION  OF  THE  "LEARN"  PROGRAM 

4 

The  LEARN  program  is  an  interpreter  for  computer-aided  instruction  scripts,  which  is 
furnished  with  the  Version  7  UNIX  system.  The  program  as  furnished  had  several  bugs  due 
either  to  transcription  or  to  changes  made  locally  to  our  system,  mostly  the  latter.  These 
errors  have  been  corrected  and  the  changes  needed  to  make  the  scripts  correspond  with  our  local 

system  have  been  completed  for  the  most-used  scripts.  A  new  script  to  familiarize  a  user  with 

1 

the  SDC  documentation  system  has  been  added. 

The  LEARN  program  interprets  scripts  designed  to  introduce  features  of  UNIX  and  to  allow 
users  to  practice  using  programs  and  commands  of  the  system  without  needing  someone  to  guide 
and  answer  questions.  The  program  is  furnished  with  six  scripts  covering  commands  and  the 
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UNIX  file  system  (two  scripts),  the  editor,  the  program  for  formatting  mathematics,  the  Bell 
Telephone  Laboratories  (BIT.)  documentation  formatting  routines,  and  the  language.  The 
script  on  the  B  IT.  documentation  has  been  rewritten  to  apply  to  the  SIX'  documentation  system. 

The  LEARN  program  presents  the  explanatory  material  from  the  script  on  the  student's 
terminal  ami  allows  the  setting  up  of  files,  special  commands,  and  the  collection  of  the  student's 
input  and/or  output.  t  he  student  is  posed  a  problem  which  may  require  some  commands  to  be 
used  or  may  require  a  specific  output  be  generated.  Other  problems  require  the  answer  to  a 
question  of  the  nature  "answer  N,"  where  N  is  a  number  dependent  on  the  question  posed  in  the 
script.  Other  script  questions  require  yes  or  no  answers.  Depending  on  the  answer,  the  score 
is  incremented  or  decremented  (between  zero  and  ten)  and  the  next  lesson  is  selected  based  on 
the  student's  current  score;  then  the  lessons  can  be  structured  in  a  tree  with  multiple  "tracks," 
depending  on  the  student's  performance. 

With  the  prospect  of  a  wider  community  of  users  of  our  system  in  the  future,  the  LEARN 
program  provides  a  valuable  capability  for  a  new  user  to  become  familiar  with  the  use  of  the 
UNIX  system  without  needing  a  tutor  to  "hand-hold"  the  new  user.  Further  additions  and  cor¬ 
rections  of  any  remaining  problems  in  the  scripts  will  be  undertaken  as  time  permits.  Com¬ 
ments  from  users  are  solicited  to  help  resolve  any  remaining  problems. 

A.  G.  Gann 
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II.  SK1SMIC  ALCORITIIMS 


A.  iI)i:\'ril'ICATI()N  OK  LOCAL  KVLNTS  IN  ACTOMATIC  SKl.VAL  DLTLC  TIO.Nb 

Ot'  tin-  17  currently  installed  SKO/ASKO  stations,  all  but  5  are  in  reninna  of  moderate  to 
high  seismic  activity.  It  is  thus  inevitable  that  a  l;i rye  proportion  of  the  signals  triggering  the 
automatic  short-period  (SI1)  event  detector  will  be  from  local  or  regional  events.  The  signals 
from  sources  at  such  distances  are  of  much  higher  frequency  than  teleseisms  but  possess  suf¬ 
ficient  energy  in  the  0.5-  to  2-llz  band  used  for  detection  to  trigger  the  SRO  detector.  In  an 
earlier  SATS*  it  was  estimated  that  over  50  percent  of  the  detections  made  at  the  Albuquerque 
SRO  appeared,  from  their  high  frequency  content,  to  result  from  events  at  distances  of  probably 
less  than  10°.  In  the  absence  of  any  complete  local  bulletin  for  the  Western  F.S.  (events  being 
located  on  a  state  basis,  with  varying  levels  of  effort),  this  could  not  be  verified. 

In  Japan,  an  extensive  seismic  network  is  operated  by  the  Japan  Meteorological  Agency 
(J  MA)  which  prepares  a  seismological  Bulletin  covering  Hokkaido,  Honshu,  and  bhikoku,  with 
less  complete  coverage  for  Kyushu  and  the  closer  islands  of  the  Ryukyu,  Bonin,  and  Kurile 
chains.  For  the  time  interval  10-21  August  1978,  we  have  used  both  the  PDF  and  J.MA  Bulletins 
to  attempt  to  identify  the  detections  made  at  Matsushiro  ASRO  LMA.lt)).  We  have  found  that  the 
frequency  content  of  the  detection  can  be  used  to  determine,  with  a  fair  degree  of  success,  which 
signals  are  from  events  at  distances  of  less  than  500  km.  This  ability  to  flag  detections  as  being 
from  close  sources  can  be  of  extreme  value  in  an  automatic  association  scheme:  large  signals 
can  be  used  to  provide  an  initial  location,  and  smaller  ones  can  reasonably  be  assumed  to  be 
from  events  which  will  not  be  detectable  at  other  stations  of  a  sparse  global  network. 

Figure  1 1  - 1  shows  the  distribution  with  magnitude  of  the  events  given  in  the  JMA  and  PDF 
Bulletins,  the  shaded  portions  indicating  detection  at  MAI O.  An  additional  60  events,  mostly 
in  the  Western  t'.S.  and  not  assigned  a  magnitude,  were  reported  in  the  PDF;  none  were  de¬ 
tected  at  MA.rO.  On  the  basis  of  such  a  small  sample,  little  can  he  said  about  the  magnitude 
thresholds  of  either  catalog  or  the  detection  capability  of  MAJO.  In  the  preceding  SATS,2  a 
50-percent  detection  threshold  of  5.0  was  obtained  for  MAJO.  This  performance  could  cer¬ 
tainly  be  improved  by  lowering  the  detector  threshold,  but  this  could  only  be  achieved  at  the  un¬ 
desirable  cost  of  greatly  increasing  the  number  of  local  events  detected.  The  JMA  Bulletin  is 
certainly  incomplete,  as  it  shows  no  increase  in  the  number  of  events  in  decreasing  magnitude. 

Of  the  129  detections  made  at  MAJO,  40  could  be  associated  with  PDF  events  at  distances 
of  10°  or  greater,  34  with  JMA  events  at  distances  of  less  than  10°,  and  55  were  unassociated. 
Sample  signal  spectra  of  local  and  teleseismic  events  are  shown  in  Fig.  11-2,  together  with  the 
background -noise  level  which  several  independent  samples  showed  to  be  very  stationary.  As 
expected,  signals  from  local  events  were  dominated  by  high  frequencies  (2  to  5  Hz),  while  tele¬ 
seisms  were  predominantly  lower  frequency  (0.3  to  2  Hz).  In  both  types  there  was  a  very  large 
variation  in  spectral  shapes  from  event-to-event,  which  is  not  encouraging  for  detector  design 
based  upon  spectral  shape  of  signal  and  noise.  There  was  some  indication  that  some  of  the  local 
signals  were  perilously  close  to,  if  not  actually,  aliased. 

We  have  attempted  to  discriminate  between  local  and  teleseismic  signals  on  the  basis  of 
their  frequency  content.  Both  frequency-  and  time-domain  measurements  were  attempted,  as 
was  the  dominant  period  at  peak  signal  amplitude.  The  most  successful  of  these  was  the  ratio 
of  the  powers  Z  Sj.(t)^  in  the  signal  sit)  after  bandpass  filtering  into  high-  and  low-frequency 
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TABLE  M-l 
NETWORK  III 


Station 

Code 

Latitude 

Longitude 

Noise 

Level 

Travel-Time 

Anomaly 

mb 

Bios 

Location 

of* 

-13.91 

-171.78 

40.00 

0.21 

-0. 10* 

New  Zealand 

ole 

82.48 

-62.40 

5.00 

-0.58 

-0.04 

Canada 

anmo 

34.95 

-106.46 

2.00 

0. 19 

-0.20 

United  States 

onto 

39.90 

32.78 

2.00 

0.10 

0.10* 

Turkey 

are 

-16.46 

-71.49 

7.00 

-0.23 

-0.10* 

Peru 

asp 

-23.68 

133.90 

3.00 

-0.  75 

-0.05 

Austral  io 

8DF 

-15.66 

-47.90 

3.00 

-0.36 

0.10* 

Brazil 

bng 

4.44 

18.55 

1.00 

-1.25 

-0.07 

Central  African  Republic 

boco 

4.59 

-74.04 

2.00 

1.21 

0.04* 

Columbio 

bod 

57.85 

1 14.18 

5.00 

-0.85 

0.10* 

USSR 

bul 

-20.14 

28.61 

4.00 

-0.72 

-0.07 

Rhodesia 

chfo 

18.79 

98.98 

2.00 

-0.65 

0.00* 

Thailand 

col 

64.90 

-147.79 

5.00 

-0.51 

0.01 

United  States 

com 

16.25 

-92.13 

25.00 

0.79 

0.00* 

Mexico 

dag 

76.77 

-18.77 

12.00 

-0.59 

-0.02 

Denmark 

EKA 

55.33 

-3.16 

8.00 

0.12 

0. 19 

United  Kingdom 

elf 

53.25 

86.27 

5.00 

-0.72 

0.10* 

USSR 

goco 

45.70 

-75.48 

2.00 

-0.21 

0.10* 

Co no do 

GBA 

13.60 

77.44 

15,00 

-0. 16 

0.04 

Indio 

GRF 

49.69 

11.21 

2.00 

0.  18 

0.24 

Federal  Republic  *>f  Germany 

HFS 

60.13 

13.70 

1,00 

-0.51 

0.05 

Sweden 

ifr 

33.52 

-5.13 

8.00 

0.32 

0.00* 

Morocco 

JYSA 

62.17 

24.87 

2.00 

0. 00* 

0.20* 

Finland 

khc 

49.13 

13.58 

3.00 

-0.32 

0.10 

Czechoslovakia 

kic 

6.36 

-4.74 

3.00 

-0.85 

-0.05* 

Ivory  Coast 

KSRS 

38.00 

128.00 

3.00 

0.00* 

0.10* 

South  Korea 

LAO 

45.68 

-106.22 

0.40 

-0. 1 1 

-0.10 

United  States 

moio 

36.31 

59.49 

2.00 

0.90 

0.00* 

Iran 

mot 

36.54 

138.21 

10.00 

-0.54 

0.00* 

japan 

mow 

-67.60 

62.88 

10.00 

-0.11 

0.11 

Austro)  ia 

mbc 

76.24 

-119.36 

6.00 

-0.41 

0.14 

Canada 

NAO 

60.82 

10.83 

0.80 

-0.75 

-0.09 

Norway 

nie 

49.42 

20.32 

5.00 

0.28 

-0.02 

Poland 

niko 

-1.27 

36.80 

2.00 

1.56 

-0.20* 

Kenya 

obn 

55.12 

36.57 

6.00 

-0.47 

0.20* 

USSR 

pns 

-16.27 

-68.47 

3.00 

0. 24 

-0.08 

Bol  ivio 

que 

30.18 

66.95 

10.00 

0.25 

0.10* 

Pakistan 

sbo 

-77.85 

166.76 

30.00 

0.69 

0. 10* 

New  Zealand 

shio 

25.57 

91 .88 

2.00 

-0.54 

0.11 

India 

sno 

-70.32 

-2.33 

>6.00 

-0.13 

-0.05* 

South  Africo 

spa 

-90.00 

0.00 

8.00 

0.02 

0.10* 

United  States 

sve 

56.80 

60.63 

6.00 

-0. 18 

0.20* 

USSR 

tom 

22.79 

5.52 

4.00 

-0.  19 

0.00* 

Algeria 

tato 

24.98 

121.49 

4.00 

0.00* 

0.00* 

Ta  iwan 

-30.  17 

-70.80 

12.00 

-0.30 

-0.10* 

Chile 

wel 

-41.29 

174.78 

30.00 

-0.31 

-0.10* 

New  Zeolond 

wes 

42.38 

-71.32 

15.00 

0.35 

0.10* 

United  States 

WRA 

-19.95 

134.35 

2.00 

-1.06 

0.00* 

A  ultra!  ia 

yak 

62.02 

129.72 

5.00 

-0.89 

0.10* 

USSR 

YKA 

62.49 

-114.60 

3.00 

-0.76 

-0.07 

Co no do 

. 

Esfimofed  using  neighboring  stations. 


bands.  The  apparently  equivalent  frequency-domain  analog  (spectral  ratios)  was  less  successful 
probably  due  to  windowing  and  truncation  effects. 

Optimum  distance  discrimination  was  found  to  be  realized  by  the  ratio  of  signal  power  in  the 
1-  to  8-llz  band  to  that  in  the  0.3-  to  l-!lz  band.  The  ASRO  instrument  response  peaks  at  about 
2.5  Hz  and  falls  off  rapidly  below  0.5  Hz  and  above  6  Hz,  so  this  is  very  close  to  the  ratio  of 
power  at  3  Hz  to  that  at  1  Hz.  Figure  II - 3  shows  this  ratio  as  a  function  of  distance  A.  It  can 
be  seen  that  all  associated  signals  with  a  ratio  greater  than  0.1  are  at  distances  of  5°  or  less, 
and  that  all  those  with  a  ratio  less  than  0.03  are  more  than  10°  distant  from  lUAJO.  Between 
these  two  values  there  is  a  "gray"  area. 

In  Fig.  II  —  3,  events  at  depths  of  less  than  70  km  are  designated  by  S,  and  those  at  greater 
depths  by  D.  A  presumed  explosion  in  Novaya  Zemlya  is  denoted  by  X.  The  intrinsically 
higher-frequency  nature  of  teleseisms  from  deep  events  can  be  seen.  At  closer  distances,  one 
would  expect  this  character  to  be  negated  to  some  extent  by  the  low  Q  structure  beneath  Japan; 
this  is  not  apparent.  While  the  distance  discriminant  probably  works  because  of  the  higher  at¬ 
tenuation  at  high  frequencies,  seismic  scaling  effects  cannot  be  discounted  since  the  closer 
events  are  mostly  much  smaller  in  magnitude  than  the  teleseism.  The  events  at  distances  of 
less  than  10°  which  fall  into  the  "gray"  area  are  indeed  among  the  largest  at  close  distances. 

Of  the  55  unassociated  signals,  all  but  15  had  ratios  of  0.1  or  greater,  indicating  that  they 
were  almost  certainly  from  local  events.  It  appears  that  local  events  recorded  at  MAJO  can 
be  readily  identified  on  the  basis  of  their  frequency  content,  and  that  such  identification  can 
certainly  be  applied  to  other  stations  beneath  which  Q  is  probably  higher  than  for  MAJO.  The 
broader-band  response  of  the  NSS  instrumentation  should  make  such  identification  even  easier, 
and  we  propose  to  test  this  as  soon  as  the  data  and  (more  importantly)  local  bulletins  become 


available. 


R.  G.  North 
M.W.  Shields 


B.  A  TEST  OF  AUTOMATIC  ASSOCIATION  ALGORITHMS  USING  SYNTHETIC  DATA 

Lists  of  synthetic  arrivals  at  a  global  network  of  stations  have  been  compiled,  and  used  to 
examine  some  of  the  characteristics  of  automatic  and  manual  methods  for  event  association  and 
location.  The  lists  are  generated  by  first  computing  a  set  of  synthetic  seismic  events,  realis¬ 
tically  distributed  in  space,  time,  and  magnitude,  and  then  computing  arrivals  at  a  prescribed 
global  network  of  stations  using  a  variety  of  acceptance  criteria  which  attempt  to  model  scat¬ 
tering  in  the  Earth,  regional  variations  in  travel  times,  station  detection  capability,  and  op¬ 
erator  error.  Details  of  this  procedure  have  been  given  earlier.^ 

The  network  of  stations  used  is  the  Network  III(SP)  listed  in  document  CCD/558,  issued  by 
the  Committee  on  Disarmament  (CD)  Group  of  Experts  on  9  March  1978.  The  stations  in  this 

network  are  listed  in  Table  II —1  which  also  includes  station  noise  levels  (as  listed  in  CCD/558), 

3  4 

station  travel-time  anomalies  (Dziewonski  ),  and  station  amplitude  biases.  Where  travel-time 
anomalies  or  amplitude  biases  were  not  available,  they  have  been  estimated  from  values  at 
neighboring  stations. 

Three  arrivals  lists,  each  corresponding  to  one  day,  were  generated  for  this  network,  and 
the  principal  characteristics  of  these  lists  (referred  to  as  A,  B,  and  C)  are  shown  in  Table  II-2. 
List  A  contains  only  larger  events,  while  list  C  includes  events  down  to  m^  =  2.1.  It  appears 
that  an  average  day  of  seismicity  will  lead  to  about  1000  arrivals  at  a  network  such  as  the  one 
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TABLE  II- 

CHARACTERISTICS  OF  TH 

2 

E  ARRIVALS  LISTS 

List  A 

List  B 

List  C 

Period  Covered 

1  day 

1  doy 

1  day 

Event  List  Complete  Down  to 

4.3 

2.7 

2.1 

Total  Number  of  Events  in  List 

32 

862 

3802 

Number  of  P  Arrivals  at  Network  III 

500 

840 

814 

Number  of  pP  Arrivals 

93 

89 

69 

Number  of  S  Arrivals 

37 

60 

30 

Number  of  Events  with  ot 

5 

7 

6 

Least  25  Arrivals 

Number  of  Events  with  at 

19 

26 

15 

Least  10  Arrivals 

Number  of  Events  with  at 

28 

41 

35 

Least  5  Arrivals 

Number  of  Events  with  at 

30 

59 

51 

Least  3  Arrivals 

Number  of  Events  with  at 

30 

172 

302 

Least  1  Arrival 

Number  of  Events  with  0  Arrivals 

2 

690 

3500 

Number  of  Local  Arrivals  (<5°) 

3 

58 

182 

Number  of  Non-local 

497 

782 

632 

Arrivals  (>5°) 
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used  here,  and  this  is  useful  for  sizing  the  task  of  data  analysis  for  such  a  network.  Natural 
variations  in  seismicity  may  lead  to  substantial  variations  from  this  average,  particularly  when 
earthquake  swarms  or  aftershock  sequences  are  present. 

The  three  arrivals  lists  were  each  analyzed  by  two  automatic-association  algorithms  cur¬ 
rently  in  use  in  the  U.S.  The  first,  which  we  refer  to  as  Procedure  1,  utilizes  P-wave  arrival 
times  and  depth-phase  arrival  times  (the  latter  are  not  used  for  association,  but  are  included 
in  the  computation  of  the  hypocenter  parameters).  To  define  an  event,  5  arrivals  are  required, 
and  array  data  are  not  used.  The  second,  Procedure  2,  does  not  use  depth-phase  information, 
but  does  utilize  array  estimates  of  azimuth  and  distance  for  trial  epicenters.  In  this  case,  a 
minimum  of  3  arrivals  is  required  to  define  an  event.  Event  solutions  for  3  or  4  arrivals 
are  computed  by  restraining  the  depth  to  zero. 

The  event  lists  generated  by  these  algorithms  were  compared  with  the  original  event  lists. 
The  quality  of  the  event  solution  was  defined  in  terms  of  the  accuracy  of  the  epicenter  location. 
If  6  is  the  sum  of  the  absolute  values  of  the  errors  in  latitude  and  longitude,  then  the  event 
quality  is  assigned  as  follows: 


6<  0.5“ 

Good 

0.5“  <64 

1.0“ 

Fair 

1.0“  <64 

10“ 

Poor 

6  >  10“ 

Bogus  event 

The  results  from  Procedure  1  are  summarized  in  Fig.  II-4.  Detection  and  event  quality 
are  excellent  for  m^  ;>.  4.5.  Below  this  magnitude,  poor  solutions  and  missing  events  begin  to 
dominate.  In  this  case,  5  bogus  events  were  generated  by  the  algorithm,  each  with  5  arrivals. 

Procedure  2  gave  improved  results  (see  Fig.  II-5).  Wore  events  were  detected  for  m^  <  4.0, 
and  the  quality  of  solutions  in  the  range  4.0  to  4.5  is  improved.  A  substantial  number  of  miss¬ 
ing  events  is  still  present,  since  the  original  lists  contained  many  events  with  3  or  4  arrivals. 
Using  this  algorithm,  13  bogus  events  were  generated  from  groups  of  3  or  4  arrivals,  and  1 
bogus  event  was  formed  from  5  arrivals. 

In  order  to  compare  these  automatic  procedures  with  the  capability  of  a  human  analyst,  an 
experienced  analyst  was  asked  to  carry  out  manual  association  on  the  same  data.  At  least 
3  arrivals  were  required  to  define  an  event.  The  results  are  shown  in  Fig.  II -6 .  The  human 
analyst  performance  is  substantially  better  than  the  automatic  procedures.  Only  a  small  num¬ 
ber  of  events  was  missing,  though  misassociations  did  lead  to  a  number  of  poor  solutions.  Three 
bogus  events  were  obtained  -  two  with  3  arrivals  and  one  with  5  arrivals. 

Several  conclusions  can  be  made  from  this  study.  Automatic  procedures  are  still  quite 
inferior  compared  with  a  human  analyst.  Additional  research  is  needed  in  an  attempt  to  im¬ 
prove  the  performance  of  the  automatic  algorithms.  Also,  the  choice  of  criterion  for  event 
definition  has  a  large  impact  on  the  final  event  list.  Network  detection  capability  is  substan¬ 
tially  improved  when  events  with  3  or  4  arrivals  are  included,  but  this  improvement  is  ac¬ 
companied  by  an  increase  in  the  number  of  missing  and  bogus  events.  Finally,  even  an  ex¬ 
perienced  analyst  is  not  perfect  and,  in  practice,  there  will  always  be  some  (hopefully  few) 
bogus  events  in  the  final  event  list. 


ftl.  A.  Chinnery 
M.  A.  Tiberio 


A  COMPARISON  OF  RELATIVE  AND  ABSOLUTE  EARTHQl  AKE 
I.OCATIONS  IN  AN  ISLAND  ARC  REGION 


Quantifying  the  effects  of  laterally  heterogeneous  Earth  structure  should  be  easiest  in 
regions  of  active  subduction  where  there  are  many  geophysical,  geochemical,  and  geological 
studies  that  claim  to  have  observed  such  effects.  Here  we  compare  earthquake  epicenters  from 
the  Bulletin  of  the  international  Seismological  Center  (ISC)  with  relative  epicenters  computed 
from  the  same  arrival-time  data.  The  relative  epicenters  are  more  clustered  and  show  an  ap¬ 
parent  westward  offset  of  about  12  km. 

The  relative  location  scheme  has  been  described  i  arlier.  ’ ^  In  the  previous  applications  of 
this  scheme,  relative  depth  was  of  as  much,  if  not  more,  interest  than  the  relative  epicenters 
because  any  accurate  measure  of  source  depth,  whether  it  be  relative  or  absolute  depth,  has  a 
potential  application  to  the  discrimination  problem.  However,  in  island  arc  regions  and  other 
regions  that  are  predominantly  water  covered,  an  accurate  epicenter  may  be  sufficient  to  identify 
an  event  as  an  earthquake,  rather  than  an  underground  explosion,  if  its  epicenter  is  in  deep  water. 

Accurate  epicenter  determinations  in  island  arc  regions  require  techniques  that,  to  some 
extent,  are  transparent  to  strong  lateral  variations  in  crust  and  upper-mantle  structure  exem¬ 
plified  by  the  inclined  zones  of  mantle  seismicity.  The  master  earthquake  technique  used  here 
is  expected  to  be  less  sensitive  to  lateral  heterogeneity  than  the  standard  bulletin  location  tech¬ 
nique,  because  the  location  problem  is  reduced  to  a  near-source  problem.  Here,  the  near¬ 
source  region  is  defined  by  a  hemisphere  of  about  100  km  radius  centered  on  the  master  earth¬ 
quake.  To  eliminate  trade-offs  between  relative  depth  and  relative  origin  time  and  epicenter, 
the  activity  is  assumed  to  occur  at  the  depth  of  the  master. 

Shallow  activity  near  the  southern  terminus  of  the  Philippine  Trench  is  well  suited  for  this 
comparison  because  of  the  intensity  of  shallow  activity  and  the  structural  complexity  which  might 
impart  systematic  errors  to  ISC  locations.  About  200  earthquakes,  of  which  only  5  were  reported 
with  a  mantle  depth,  were  located  relative  to  the  major-size  earthquake  (Mg  7.8)  of  2  December 
1972.  The  source  region  is  shown  in  Fig.  II-7.  The  ISC  a.  '  relative  epicenters  are  shown  in 
Figs.  II — 8  and  II-9,  respectively.  Clearly  the  relative  epicenters  are  more  clustered  and,  in 
particular,  the  clustering  reveals  at  least  three  aseismic  regions  surrounded  by  comparatively 
intense  activity  (regions  A,  B,  and  C  in  Fig.  II-9).  In  addition  to  the  clustering,  there  is  an 
apparent  shift  of  the  relative  epicenters  to  the  west  by  about  12  km  relative  to  the  ISC  epicenters. 
If  this  is  interpreted  as  a  systematic  error  in  the  ISC  locations,  then  the  ISC  location  for  the 
master  must  be  taken  as  correct.  We  cannot  prove  such  an  assertion  but  merely  point  out  that 
the  azimuthal  coverage  for  the  epicenter  determination  of  the  master  is  as  good  as  can  be  ex¬ 
pected  for  this  region  and,  consequently,  an  accurate  determination  is  possible  in  spite  of  the 
structure  complexity  of  the  region  because  the  epicenter  is  determined  more  on  geometric  con¬ 
siderations  than  on  structural  ones. 

T.  J.  Fitch 
R.  Catchings 

D.  SOME  SUGGESTIONS  FOR  IMPROVING  THE  ACCURACY 

OF  HYPOCENTER  DEPTH  ESTIMATES 

The  purpose  of  this  report  is  to  describe  some  problems  encountered  in  the  use  of  pP  and 
sP  for  depth  determination  of  shallow  focus  oceanic  earthquakes.  Some  ways  of  overcoming 
these  problems  are  suggested. 


TABLE  11-3 

DATES,  LOCATIONS,  DEPTHS,  AND  MAGNITUDES  OF  EVENTS  ANALYZED 


Time 

Longitude 

H 

mb 

M 

s 

1 0/23/64 

01:56:05. I 

19.80 

56.11 

25 

43 

31 

6.2 

ED 

08/19/77 

05:08:41.6 

-11.15 

118.39 

18 

54 

33 

6.  1 

C8/20/77 

19:16:32.7 

119. 13 

17 

33 

33 

6.0 

6.1 

08/25/77 

18:05:10.8 

-10.74 

119.26 

1 1 

51 

33 

6.  1 

6.0 

08/26/77 

08:26:37.5 

-10.66 

119.30 

K 

33 

33 

5.6 

09/02/77 

10:36:28.3 

-11.04 

119.14 

74 

33 

6.0 

09/23/77 

05:57:55.6 

-11.20 

118.21 

34 

33 

6.0 

1 0/07/77 

12:10:43.7 

-09.96 

1)7.31 

29 

33 

5.9 

10/16/77 

21:09:17.7 

-09.73 

117.11 

39 

33 

5.6 

5.8 

SPZ  and  LPZ  WWSSN  seismograms  of  one  Atlantic  and  seven  Western  Pacific  earthquakes 
provide  the  database  for  the  analysis.  Table  II-3  lists  source  parameters,  and  compares  depths 
found  in  this  study  with  those  published  in  the  TDK  and  ISC  liulletins.  1  he  analysis,  following  a 
procedure  described  by  Forsyth/’  consists  of  a  search  for  phase  arrivals  within  10  s  after  the 
P  arrival  on  each  in  a  suite  of  SPZ  seismograms.  If  the  earthquake  is  of  sufficient  magnitude 
(M  >  5.  5),  the  P  waveforms  on  the  I.PZ  seismograms  are  also  examined. 

SPZ.  seismograms  from  the  earthquake  of  20  August  1977  are  shown  in  Fig.  11-10.  In  each 
case,  the  important  phases  are  marked  on  the  SPZ  seismogram  of  the  station  named  on  the  left. 
The  seismograms  are  lined  up  so  that  the  zero  line  represents  the  first  arrival  on  each  seis¬ 
mogram.  In  the  simplest  case,  one  would  expect  to  see  the  predominant  phases:  the  direct  P 
arrival,  the  P  wave  reflected  from  the  rock-water  interface  just  above  the  hypocenter  (pP), 
and  the  S  wave  reflected  as  P  from  the  same  area  (sP)  plus  later  phase  from  reverberations 
in  the  water  layer.  Furthermore,  the  pP-P  and  sP-P  time  differences  are  functions  mainly  of 
depth;  at  crustal  depths,  these  times  vary  only  a  few  tenths  of  a  second  over  the  distance  range 
from  30°  to  90°.  Therefore,  these  phases  would  nearly  line  up  on  a  diagram  such  as  Fig.  11-10. 

Figure  1 1 -1 0  departs  from  the  ideal  and  thus  illustrates  several  problems  in  identifying 
depth  phases.  First  of  all,  an  earthquake  is  sometimes  closely  preceded  by  a  much  smaller 
event,  called  a  precursor.  .Second,  pP  and  sP  are  often  not  particularly  prominent  on  SPZ 
seismograms.  Third,  water  reflections  are  often  mistaken  for  pP  and  sP.  Finally,  the  ar¬ 
rival  of  a  later  phase  may  be  lost  in  the  coda  of  an  earlier  one. 

The  precursor  problem  is  common  enough  so  that  it  must  always  be  anticipated.  Figure  11-10 
shows  a  precursor  on  several  seismograms.  It  is  recognized  as  such  because,  over  the  entire 
110°  range  of  azimuths,  its  amplitude  remains  much  smaller  than  that  of  the  next  succeeding 
arrival  (the  direct  P  arrival  of  the  main  shock).  Carefully  monitoring  the  relative  amplitudes 
of  the  early  phases  should  be  an  aid  in  identifying  these  small  events.  For  example,  Kanamori 
and  Stewart7  judge  the  Gibbs  Fracture  Zone  earthquake  of  16  October  1074  to  be  an  11-km-deep 
event  with  a  precursor.  They  identify  the  precursor  on  the  basis  of  low-amplitude  first  arrivals 


TABLE  11-4 


EXPECTED  X-P  TIMES  FOR  FOUR  PHASES  ASSUMING  MODELS 
OF  TABLE  11-5  AND  VERTICAL  RAY  PATHS 
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1 

pWwP  | 

! 

swwP 

Earthquake  Depth 
(Belcw  Sea  Floor) 

(km) 

10.8 
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16.8 

10 
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15 

14.2 
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20 
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20. 1 

10 

15.3 
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15 
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TABLE  11-5 

PARAMETERS  OF  OCEANIC  BASIN  AND  TRENCH  MODELS* 


Thickness 

(km) 


Trench 


Water 

Sediment  Layer  1 
Sediment  Layer  2 
Crust 
Mantle 


*  Table  3  from  Ward? 


P 

Velocity 

(km/s) 


Velocity 

(km/s) 


Density 

(g/cm3) 


IP*!#*  jm‘- 


on  all  seismograms.  An  11 -km  focal  depth  is  reasonable  for  an  earthquake  located  on  a  trans¬ 
form  fault.  The  P  wave  of  the  main  shock  follows  the  precursor  by  11  to  12  s  and  is  identified 
in  the  ISC  Bulletin  as  pP.  Depth  is  estimated  to  be  41  km,  unreasonably  deep  for  an  earthquake 
on  a  transform  fault. 

The  last  three  problems  mentioned  above  are  worst  for  earthquakes  with  m^  ■  5.5.  The 
lack  of  prominence  of  pP  and  sP  on  SPZ  seismograms  arises  from  two  causes.  First,  short- 
period  energy  is  sensitive  to  finer  detail  of  crustal  structure  resulting  in  a  greater  number  of 
reflections.  If  m^  >  5.5,  more  energy  is  available  for  transmission  through  and  reflection  from 
various  crustal  layers,  increasing  the  number  of  phases  recorded  on  the  seismogram.  Further¬ 
more,  if  a  thick  layer  of  soft  sediment  lies  above  the  harder  crust,  the  reflection  coefficients 
for  pP  and  sP  are  much  reduced  relative  to  reflections  from  other  layers.  The  second  cause 
is  variation  of  amplitude  due  to  radiation  pattern.  The  MAT  SPZ  seismogram  of  Fig.  11-10 
illustrates  this  best.  The  focal  mechanism  deduced  for  this  earthquake  predicts  a  pP  relative 
amplitude  of  nearly  0.0  at  MAT,  a  prediction  justified  by  the  insignificant  pP  phase  on  the 
seismogram. 

Reflected  phases  most  often  confused  with  pP  and  sP  for  shallow  oceanic  events  are  re- 

g 

flections  from  the  water-air  interface  called  pwP  and  swP  (Chen  and  Forsyth  ).  These  phases 
are  quite  prominent  on  short-period  seismograms,  and  it  was  found  in  Forsyth's  and  in  this 
study  that  their  misidentification  often  leads  to  overestimation  of  hypocentral  depth  by  factors 
of  2  or  3  in  bulletins.  Table  IX  —4  presents  expected  phase-P  times  for  four  water-air  reflected 

9 

phases  based  on  the  oceanic  basin  and  trench  models  of  Table  11-5  (Ward  ). 

The  water-reflection  problem  combines  with  the  coda  problem  as  the  major  causes  of  in¬ 
correct  focal -depth  estimates.  On  seismograms  of  large,  shallow  oceanic  shocks,  the  pP 
arrival  is  often  obscured  by  the  coda  of  P.  For  a  10-km-deep  quake,  the  pP-P  time  is  only 
3  s  and  sP-P  is  4.3  s  on  teleseismic  records.  If  the  quake  magnitude  is  large,  all  three  phases 
and  their  codas  are  inseparable  and  thus  are  reported  as  one.  The  next  identifiable  phase  at 
10  s  after  P  is  pwP.  Taking  this  to  be  pP  results  in  a  depth  estimate  of  about  33  km,  a  three¬ 
fold  overestimate.  Nearly  all  the  depth  discrepancies  between  this  study  and  the  ISC  Bulletin 
can  be  traced  to  misidentification  of  pwP  or  swP  as  pP. 

The  chaotic  nature  of  the  SPZ  seismogram  for  a  large  shallow  oceanic  earthquake  makes 
it  important  that  the  early  part  of  the  LPZ  WWSSN  seismogram  be  analyzed  when  available. 

The  impulse  response  of  the  LPZ  instrument  is  shown  in  Fig.  11-11.  The  analysis  assumes 
that  for  shallow  events  the  LPZ  P  waveform  is  the  convolution  of  impulsive  P,  pP,  and  sP 
arrivals.  Phases  pP  and  sP  can  often  be  identified  as  points  where  the  waveform  departs  from 
the  impulse  response  to  a  single  phase  arrival.  These  two  phases  are  relatively  more  prom¬ 
inent  on  the  LPZ  than  the  SPZ  record  because  pwP  and  swP  are  mainly  short-period  phases 
(Forsyth^).  On  LPZ  seismograms  for  relatively  deep  events  (e.g..  Fig.  11-12,  depth  14  to 
16  km),  a  pP  or  sP  arrival  will  be  seen  as  a  shelf-like  offset  of  the  early  part  of  the  waveform. 
On  LPZ  seismograms  from  a  very  shallow  event  (e.g..  Fig,  11-13,  depth  <  10  km),  the  P  wave¬ 
form  will  be  a  smooth  convolution  of  three  impulses  and  will  not  resemble  the  LPZ  WWSSN 
impulse  response.  The  deduction  follows  that  SPZ  arrivals  with  phase-P  times  of  10  s  or 
longer  are  probably  water  reflections. 

For  earthquakes  with  magnitudes  less  than  about  =  5.5,  the  P  wave  signals  are  too 
weak  to  be  seen  on  most  LPZ  records  but  the  SPZ  record  is  usually  easier  to  analyze.  Incor¬ 
rect  focal-depth  estimates  are  most  likely  to  be  the  result  of  failure  to  recognize  a  precursor 
or  a  low -amplitude  pP. 


Two  other  phenomena  .should  also  ho  noted  as  qualitative  indicators  of  shallowness,  first, 
on  the  1  IV.  record,  a  signal  pictured  m  fig.  II-l  1  and  called  "ringing  P"  by  Ward  '  is  a  prob¬ 
able  indicator  that  most  of  the  fault  break  is  in  the  crust.  I. ike  Ward,  the  amplitude  of  the  ring¬ 
ing  I’  waie  we  observed  did  not  vary  systematically  with  focal  depth.  Second,  to  12-s  period 
surface  waves  on  si’/,  records  of  earthquakes  with  moderate  magnitudes  indicate  shallowness. 

M.  W.  Shields 

!•:,  \i:w  Ai’i’i  . it  a  i  m  ).\s  st  if  i  w  Am: 

Pour  pieces  of  software  which  recentlv  have  been  released  are  described  below. 

(It  a  power  detector.  l*resentl\,  w e  have  no  need  for  a  delet  tor;  however,  m  the  future 
we  expert  that  continuous  si'  data  will  be  fed  to  the  data  (enter  and  so  we  have  begun  testing  a 
detector.  Since  real-time  data  arc  not  available,  we  have  used  a  tape  of  continuous  Ipredetm  tor) 
SKO  output  which  has  enabled  us  to  optimize  our  detector,  the  limitations  being  that  our  data 
tiles  are  ot  tinitr  length,  lit  uses  a  sliding  window  l.k  s  long  over  wincii  the  power  is  calculated, 
com  pa  ml  with  'he  a  vr  rage  ponvi  st  .  1 1  s'  n  s  foi  noise  i  ik.n  If  s  pi  m  r  'o  tin-  t  ■■  s  I  point ,  a  s  its 
basic  detection  criterion. 

When  the  test  power  exceeds  5m  times  the  standard  deviation  of  the  average  noise,  a  de¬ 
tection  is  declared  and  the  average  noi»  •  statistic  s  are  fro/en.  In  the  same  fashion  as  the  SKO 
detector,  our  program  places  a  marker  "STAKT"  IS  s  before  the  arrival.  The  arrival  time  is 
then  marked  by  "AK\  f."  The  dolrctmr  is  declared  over  when  the  power  is  below  the  threshold 
lor  at  least  lb  S ,  at  which  pun;t  a  marker  "I.Ni)"  is  written.  I  ills  leaves  the  important  parts 
of  the  continuous  record  labeled  for  extraction  latm  .  \  bulletin  of  arrivals  and  detector  on  and 

off  times  is  printed. 

ppick:  a  specialized  version  of  dt.  Slit)  records  presently  are  received  with  little  indica¬ 
tion  ot  where  the  arrival  is  in  the  record,  for  \  vese  •  refer  or  inalvst  to  make  the  picks  e  .n 
sometimes  he  a  long  and  tedious  task.  1’piek  is  written  to  he  fast;  it  will  search  the  first  s 
of  the  SKO  record.  Ppick  gets  its  noise  averages  from  the  first  14  s  of  the  record  and  doe.- 
not  update  with  later  values.  Thus,  if  there  is  less  than  14  s  of  noise  before  the  arrival,  chances 
are  the  arrival  will  not  he  picked.  \n  example  of  ppick 's  performance  is  given  in  fig.  11-14 
which  shows  the  P  wave  arrivals  from  some  selected  Pastern  Kazakh  "events"  as  recorded  by 
the  Albuquerque,  New  Mexico  SKO.  The  markers  without  labels  are  produced  by  ppick.  Ppick 
makes  no  attempt  to  identify  the  phase.  Those  markers  labeled  ".I -15"  are  produced  bv  another 
piece  of  software  we  have.  Their  position  is  calculated  from  th  ;  station  and  event  location  as 
well  as  .1-15  travel-time  tallies.  The  markers  labeled  "11"  were  placed  there  by  myself.  The 
superior  quality  of  the  detector  picks  over  the  .1  —  15  calculated  picks  is  clearly  evident.  The  de¬ 
tector  has  the  added  bonus  that  it  requires  no  a  priori  knowledge  of  the  event. 

grmers:  a  graphics  program.  Clrmcrs  can  display  the  phase  spectrum  of  the  crosscot  - 
relation  of  two  gramfiles.  It  can  also  display  the  amplitude  spectrum  of  the  autocorrelation  of 
a  gramfile.  This  program  was  developed  at  the  request  of  one  of  our  researchers  who  was 
looking  for  a  way  to  view  the  source  generated  phase  differences  between  path  equalized  seis¬ 
mograms.  The  program,  which  works  on  I  P  and  SP  data,  was  later  expanded  to  its  present 
state. 

polypack-,  an  interactive  polynomial  calculator.  Polvpack  was  developed  during  an  instru¬ 
ment  response  study.  Essentially,  the  user  can  input  u p  to  <50  polynomials  of  up  to  3Sth  degree. 


Basic  mathematical  manipulations  such  as  addition,  subtraction,  multiplication,  division,  square 
root,  root  finding,  and  evaluation  are  available.  All  math  is  done  in  double  precision. 


M.  A.  Tiberio 
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Fig.  [1-3.  Ratio  of  signal  power  in  3-  to  8-Hz  and  0.3-  to  1-Hz  bands 
as  a  function  of  epicentral  distance  A  in  degrees. 
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Fig.  11-4.  Event  quality  for  automatic  association  as  a  function  of  true  event 
magnitude  (Procedure  I).  Two  events  with  >  6.0  are  not  shown.  Counts 
are  summed  over  lists  A,  13,  and  C. 


Fig.  II-5.  Event  quality  for  automatic  association  as  a  function  of  true 
event  magnitude  (Procedure  2).  Two  events  with  >  6.0  are  not  shown. 
Counts  are  summed  over  lists  A,  B,  and  C. 
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Fig.  II-10.  SPZ  WWSSN  seismograms  for  earthquake  of  20  August  1977, 


Fig.  II  —  1 1 .  Impulse  response  for  WWSSN  LP 
instrument  broadened  by  a  t*  of  0.5  s.  Rise 
and  fall  times  of  impulse  =  1.00  s;  width  = 
0.00  s. 
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Fig.  11-13.  LPZ  WWSSN  seismograms  for  earthquake  of  2  September  1977. 
P-waveforms  are  not  offset  like  those  of  20  August  1977,  but  deviate  from 
expected  response  to  a  single  impulse.  Earthquake  depth  <10  km. 


Fig.  11-14.  Comparison  of  Human  (H)  and 
Detector  (no  label)  made  picks.  J-B  table 
calculated  onset  times  are  also  shown  (J-B). 
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III.  GENERAL  SEISMOLOGY 


A.  S!  REACH  WAVES  GENERATED  MY  COMBINATIONS  OF  EXPLOSIVE 

AM)  FACETING  SOCRCES 

In  the  preceding  SATS,1  the  results  of  a  study  of  SRO  data  from  several  presumed  explosions 
in  Eastern  Kazakh  were  presented.  For  one  of  these  events  (7  July  1979),  we  found  that  the 
Rayleigh  waves  were  reversed  in  phase  over  the  entire  usable  period  band  (15  to  50  s)  at  all 
azimuths  and  that,  in  addition,  the  Rayleigh  waves  from  this  event  displayed  significant  azi¬ 
muthal  asymmetry  compared  with  all  the  other  events  from  the  same  small  source  region. 

Here,  we  consider  the  character  of  the  far-field  surface  radiation  generated  by  the  addition  of 
a  faulting  component  to  the  symmetric  explosion  source.  The  most  convenient  representation 
of  such  a  combination  is  given  by  the  seismic  moment  tensor,  to  which  the  radiation  is  linearly 
related. 

The  far-field  Love  (FR)  and  Rayleigh  (FR)  radiation  from  a  source  of  given  moment  tensor 
is  given  by^ 

F  |  =  { A  j  [  Mxy  cos  2  0  +  0. 5(  Myy  -  Mxx)  sin  20  ) 

+  iB^  (-Myz  cos  0  +  Mxz  sin0]}  e*1^4 

Fr  =  {AR  [Mxy  sin20  -  0. 5(Myy  -  Mxx)  cos  20  1 

+  CrMzz  4  0.5AR(Myy  +  Mxx) 

-  illR  [  Myz  sin  0  +  Mxz  cos  01}  e'”^4 

where  AR,  Bj  ,  AR,  B R,  and  CR  are  related  to  the  displacements  and  stresses  at  the  source 
depth  h.  For  very  shallow  sources  such  that  (h/X)  «  1,  the  stress  components,  entering  into 
and  Br,  are  very  small  at  a  given  wavelength  X  and  vanish  at  the  surface.  Additionally, 

CR  =  -l/3AR  close  to  the  free  surface.  Consequently,  at  shallow  depths  the  imaginary  parts 
of  Fr  and  FR  are  very  poorly  excited,  and  the  above  expressions  reduce  to 

F  ^  =ARe  ^'r^4[  Mxy  cos  20  +  0.5  (Myy  -  Mxx)  sin  20] 

and 

Fr  =  Ar  ei,T^4  [Mxy  sin20  -  0.5(Myy  -  Mxx)  cos  20 

-  l/3  Mzz  +  0.5(Myy  +  Mxx)] 

Note  that  the  source  phase,  excluding  the  e±i7r^4  terms,  can  only  be  0  or  n  since  the  imagi¬ 
nary  part  essentially  vanishes,  and  that  the  shape  of  the  spectrum  is  given  simply  by  either  A^ 
or  Ar,  regardless  of  the  azimuth  of  observation  0. 

For  an  explosive  source,  the  moment  tensor  is  given  by  equal  diagonal  components 

Mxx  =  Myz  =  Mzz  =  M  ® 
producing  azimuthally  invariant  radiation. 
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Following  Mendiguren,  the  tensor  components  for  a  fault  aligned  in  the  x-direction  of  a 
C  artesian  (x,  v,  z)  system  are  given  by: 

F* 

Vertical  strike-slip  Mxy  =  M  ,  others  zero 

Fure  dip-slip  Mxy  =  Mxx  =  0 

P 

Myy  =  -M  sin  2d 

p 

Mzz  =  +M  sin  2d 
o 

p 

Myz  =  — M  cos  2d 

F 

Mxz  =  +M  cos  2d 
o 

where  d  is  the  dip  of  the  fault  plane,  and  0  <  d  <  90  represents  thrust  faulting  and  90  <  d  <  180 
normal  faulting. 

In  the  no  tat.  on  of  Mendiguren,2  the  far-field  Rayleigh  Friuj)  and  Love  F,  (tv)  are  given  as 
follows: 

For  an  explosive  source, 

Fr  =  2/3ARMQe 


For  the  vertical  strik  '-slip  source, 

F 

Fr  =  ArM  sin  2 0 

fl  =  alm<T  cos20 

For  a  pure  dip-slip  source, 

Fr  =  l/2ARM^  sin  2d  (cos  20  -  5/3) 

F^  =  -l/2A^M^  sin2d  sin20 

where  we  have  everywhere  dropped  the  e±i,r^4  terms.  AL  is  always  positive,  and  for  shallow 
depths  Ar  is  negative  over  the  period  range  of  interest.  The  source  phase  of  the  Rayleigh 
waves  from  an  explosive  source  is  thus  n,  and  that  for  dip-slip  faulting  is  0  for  sin  2d  >  0 
(thrust)  and  n  for  sin  2d  <  0  (normal  faulting). 

We  are  particularly  interested  in  combinations  of  explosive  and  faulting  sources  which  will 
generate  phase  reversals  (i.e.,  source  phase  of  0)  at  all  azimuths  0,  compared  with  the  pure 
explosive  source.  The  strike-slip  source  generates  the  familiar  four-lobed  pattern  for  which 
the  source  phase  is  alternately  zero  or  li  and  the  source  phase  for  a  normal  fault  is  every¬ 
where  n.  Only  the  thrust  fault  has  a  source  phase  of  zero  at  all  azimuths  0.  A  combination 
of  the  explosive  and  thrust-faulting  sources  will  thus  produce  tne  desired  pnase  reversal  if 
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p 

M  is  sufficiently  large.  The  radiation  from  the  thrust  fault  is  maximized  for  a  dip  of  45*,  and 
for  this  case  the  spectrum  of  such  a  combination  is  given  by 


Fd  =  ARM®  [2/3  +  1/2Rfr(cos20  -5/3)1 


FB' 


where 


M 


FB 


M 


Love  waves  are  excited  only  by  the  thrust  faulting  and  are  given  by 
FL=TMoRFBsi"2^  • 


The  Rayleigh-  and  Love-wave  radiation  patterns  are  given  in  Fig.  III-l  for  various  values 
of  Rpg.  The  Love-wave  radiation  pattern  grows  in  amplitude  linearly  as  Rpg  increases.  The 
Rayleigh-wave  radiation  patterns  exhibit  the  interesting  behavior  of  changing  from  circular 
(explosion)  to  two-lobed  (normal  faulting)  through  four-lobed  (vertical  strike-slip  faulting)  to 
two-lobed  (thrust  faulting)  as  RpB  increases.  On  the  basis  of  the  Rayleigh-wave  radiation 
alone,  these  patterns  for  different  types  of  faulting  are  indistinguishable  from  a  combination 
of  explosive  and  thrust-faulting  sources.  Note  that  the  phase  reversal  at  all  azimuths  does 
not  occur  until  RRR  exceeds  2.  The  relative  amplitudes  of  Love  and  Rayleigh  are  a  function 
of  frequency,  and  for  a  shield  modelR  the  Love-wave  amplitudes  shown  should  be  multiplied  by 
(1.36,  1.79,  2.17,  1.95,  and  1.88)  at  periods  of  50,  40,  30,  20  and  16  s,  respectively,  to  cor¬ 
rectly  reflect  relative  Love  and  Rayleigh  amplitudes  at  these  periods. 

It  has  been  pointed  out3  that  faulting  involved  with  the  explosive  process  need  not  involve 

4 

tectonics,  and  that  fracturing  along  dipping  planes  has  been  observed  at  the  NTS.  Such  fractures 
must,  however,  have  a  preferred  dir  -ction  to  produce  the  features  of  tne  surface-wave  radiation 
observed  for  some  Eastern  Kazakh  explosions.1  Thrust  faulting  on  a  conical  surface  above  the 
explosive  source  may  be  modeled  by  considering  it  as  a  sum  (or  integral)  over  all  possible  fault 
orientations  <t>.  In  this  case,  for  a  pure  dip-slip  fault  at  dip  angle  d,  orientation  0,  we  have 

Mxy  =  l/2  sin  0  sin  2d 

Mxx  =  -sin^  0  sin  2d 
7 

Myy  =  -cos  0  sin  2d 

Mzz  =  sin  2d 


and 


Myz  =  -cos  0  cos  2d 
Mxz  =  sin  0  cos  2d 


=  0 
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The  net  result  is  a  compensated  linear  vector  dipole  (CLVD)  with  zero  Rayleigh-wave  phase 
which  does  not,  however,  generate  Love  waves  or  any  azimuthal  variation  in  the  Rayleigh-wave 
radiation.  These  last  two  features  require  a  preferred  orientation  of  faulting  which  may  lie  con¬ 
ditioned  by  either  pre-existing  stress  or  planes  of  weakness  in  the  source  region. 

We  plan  to  investigate  further  the  character  of  surface-wave  radiation,  with  particular 
emphasis  upon  the  Love  waves,  from  explosive  sources. 

The  effect  of  the  combined  explosion/thrust  faulting  upon  Mg  is  summarized  in  Fig.  HI -2 
which  shows  the  maximum,  minimum,  and  average  (over  azimuth)  for  field  Rayleigh-wave 
amplitudes  as  a  function  of  Ilpy.  The  last  may  be  considered  to  reflect  the  average  Mg  com¬ 
puted  from  a  number  of  stations,  and  can  be  seen  to  possess  a  minimum  of  one-third  that  from 
the  pure  explosion  at  R ^  B  =  °-6‘  This  corresponds  to  a  reduction  in  Mg  of  about  0.5  unit,  which 
is  certainly  significant  in  Mg-yield  calculations. 

R.  G.  North 
T.  J.  Fitch 

B.  MEASUREMENT  OF  DISPERSION  AND  ATTENUATION  OF  MANTLE  WAVES 

THROUGH  LINEARIZED  INVERSION  OF  WAVEFORM  DATA:  I.  THEORY 

AND  A  TEST  ON  SYNTHETIC  DATA 

Much  of  our  current  knowledge  on  the  distribution  of  shear  velocity  and  attenuation  in  the 
upper  mantle  comes  from  the  studies  of  surface  waves.  So-called  mantle  waves  have  periods 
greater  than  100  s  and  are  attenuated  slowly  enough  that  it  is  possible  to  observe  at  a  given 
station  several  of  their  consecutive  passages  around  the  world.  These  LP  waves  are  now  rou¬ 
tinely  detected  by  modern  instruments  for  events  of  a  surface-wave  magnitude  as  low  as  6.5. 

In  addition  to  their  important  role  in  providing  constraints  on  the  mantle  structure  in  a  depth 
range  from  100  to  400  km,  they  are  also  used  in  source  mechanism  studies.  Precise  knowl¬ 
edge  of  phase  velocities  is  needed  in  inversion  of  waveform  data  for  the  seismic  moment  tensor. 

The  principle  of  measuring  the  phase  velocity  and  attenuation  from  observations  of  mantle 
waves  is  very  simple.  The  spectral  ratio  of  two  wave  groups  whose  paths  differ  by  one  cir¬ 
cumference  yields  a  spectrum  of  transfer  function  T(w)  that  in  the  time  domain  corresponds  to 
the  shape  of  a  6-function  signal  after  one  complete  passage  around  the  Earth5 

T(w)  =  Rnf2(0))/Rn(a>)  .  (III-l) 

The  spectrum  of  T(w)  must  be  of  the  form: 

T(a0  =  exp |  2 rra  {i  [k(w)  I  -  (IH-2) 

where  k(w)  is  the  wave  number,  u(w)  is  the  group  velocity,  and  q(u>)  is  the  inverse  of  the  quality 
factor  Q.  Thus,  the  analysis  of  the  phase  of  the  T(w)  allows  the  wave  number  to  be  determined 
as  a  function  of  frequency  and  its  amplitude  yields  Q(w),  if  the  group  velocity  is  known.  The 
frequencies  of  zero  crossings  of  the  phase  spectrum  correspond  to  the  frequencies  of  free  oscil¬ 
lations  of  the  Earth. 

In  practice,  direct  application  of  Eq.  ( 1 1 1  - 1 )  leads  to  serious  difficulties,  as  the  process  of 
deconvolution,  implied  by  the  ratio  of  two  spectra  contaminated  by  the  noise  and  energy  of  over¬ 
tones,  leads  to  unstable  results  even  lor  a  limited  range  of  frequencies. 

Nakanishi  derived  the  T(w)  through  a  design  of  the  Wiener  filter,  but  it  appears  that  this 
procedure  may  also  lead  to  physically  unacceptable  results;  the  T(uj)  shown  in  his  Fig.  2a  cannot 
be  a  functional  of  a  realistic  Earth  model. 
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Wo  intend  to  make  use  of  the  fact  that  the  Ttu)  must  bo  consistent  with  a  reasonable 
I'.arth  structure,  this  is  equivalent  to  using  tlio  known  proportion  of  tho  liarth  as  a  smoothing 
filter.  Tin  most  iliroot  way  to  assuro  that  tiiis  condition  is  satisfied  is  to  ii>  rive  tin 
TU  I  through  the  process  of  porturhation  of  structural  parameters  of  a  selected  starting  liarth 
model. 

We  shall  attempt  to  find  a  TU  )  that  satisfies  the  following  least-squares  condition: 

\  ~!  K  1+7(w)  -  K,U)  •  T(w)|2  d^'  =  min  .  (Ill-  3) 

‘"'l 

Let  the  Tu)  computed  for  the  starting  model  he  T°(u).  It  is  characterized  by  a  complex 
function  K^(u‘),  closely  related  to  the  complex  wave  number  and  corresponding  to  the  expres¬ 
sion  within  the  braces  in  Kq.  (Ill  —2) .  Thus, 

T°(a)  =  exp  l-2iriaK°(w)  I  .  ( III  -  4) 

A  small  perturbation  in  the  structural  parameters  of  the  starting  model  will  lead  to  a  change 

0  7 

in  K  (w)  that  can  he  evaluated  using  the  linear  perturbation  theory  : 

6K°(o.)  =  \  G°lr,  u)  •  6m(r,  w)  dr  (III-5) 

"O' 

where  the  differential  kernels  G®  and  the  liarth  model  m  may  both  be  complex.  The  effect  of 
attenuation  on  velocity  dispersion'  may  be  incorporated  in  Kq.  ( II I  -  5) .  If  the  change  in  6m  is 
very  small,  then  we  may  write 


exp[-2ria<5K°(u!)|  ~  1  -  2iriaAK° U)  =  1  -  2iria  \  G 


The  integral  over  a  finite  frequency  range  in  Kq.  (III-?)  may  be  changed  to  a  sum  over  the 
discrete  elements  of  the  spectrum,  and  upon  substitution  of  Kq.  (til -61  the  least-squares  con¬ 
dition  becomes 


"  «„<“>, Tj°  +  2*‘«H|1(wl1T|0  j  Qj  ■  dmeir  j  2  -  min 


(III-7) 


The  inverse  problem  for  Amir)  posed  by  Kq.  IIII-7)  can  be  solved  either  in  the  data  or  pa¬ 
rameter  space  using  one  of  many  methods  described  in  the  literature.  It  is  important,  however, 
that  the  perturbation  estimated  in  each  iteration  be  small  enough  that  the  approximation  assumed 
in  Kq.  (IIL6)  is  valid.  1-or  this  reason,  it  will  be  generally  necessary  to  damp  the  inverse  heavily 
and  perform  several  iterations.  After  each  iteration,  the  T(w)  is  recalculated  exactly: 


T,i)(wl 


T1''"  exp[  ~2riaAK(-')| 


After  the  convergence  is  achieved,  it  is  prudent  to  recompute  function  Kiev)  and  the  differ¬ 
ential  kernels  for  the  new  model  and  repeat  the  entire  process. 

Another  difficulty  in  estimation  of  tho  TG  I  is  that  the  wavegroups  corresponding  to  the  in¬ 
dividual  orbits  cannot  always  lie  readily  isolated.  Glearlv,  when  the  receiver  is  close  to  the 
source  or  its  antipodes,  then  the  energy  of  odd-  and  even-numbered  wavegroups  will  overlap. 
However,  even  at  the  optimal  distance  of  ''0“,  a  SOO-s  wave  of  will  arrive  at  the  same  time 


(  i  »Mi\r  <»l'  K  y  I'  tir  tin iva.vin,  iht?  most  reliable  incasiireinenta  uf dispersion  and  at- 

ot  ?■  .inti**  w have  been  limited  to  the  period  ran^e  for  ISO  to  1U0  s  -  the  interval 
iii  Aim  h  both  Kavleirh  and  1  ove  navis  are  weakly  dispersed. 

1  tie  development  "f  the  technique  described  above,  which  does  not  require  estimation  of  the 
. , ;  it  ‘i  tral  rati.'  ir.t  r  « •  1 1 .  ■  >  <  >i  n  mi  n .  i  /  mu  the  difference  between  two  functions  of  time  or  frequency, 
allow..  i>  to  i»-  miniate  the  approach  to  the  analysts  of  mantle  waves. 

(  oiir-ide?  that  a  rmmi'  o^rMn  consists  of  a  series  of  arrivals  of  wave^roups:  Rj,  R^,  R^, 

K  .  .  (  |{  _  .  t  .  >r olutioii  t*f  such  a  seismogram  with  the  transfer  function  T  transforms  it  into 
a  serif."  R  r  K  Ks,  .  .  .  ,  K,  .  The  difference  between  these  two  series  is  Rj  *  R^,  a  signal  of 
finite  duration  t^,  d«*te t  in ined  hv  the  slowest  group  velocity  of  the  wavegroup  R-.  The  remainder 
of  the  trai  e  i identieally  zero,  if  the  transfer  function  T  cor  responds  to  that  for  the  real  Karth. 
In  practn  e,  we  shall  he  dealing  with  seismograms  of  finite  duration  trunc  ated,  say,  at  a  time  t^. 
Thu.",  the  operation  described  above  is  equivalent  to 

Mtl  -  St t)  >  T(t)  0  for  t.  ^  t<t2  .  <111-9) 

In  the  frequency  domain,  the  least-squares  condition  thus  w i XI  be 

2 

F  (S(t)l.  -  F  l-S(t)  *  T(t)l  da  =  min  (III-10) 

ll,2  t,  2 

where  F  is  the  Fourier  transform  and  the  subscript  indicates  that  the  transform  is  evaluated 
for  the  time  interval  from  tj  to  t2<  Further  procedure  is  identical  to  that  described  above  for 
a  pair  of  the  wavegroups. 

What  has  been  gained  is  that  we  can  use  the  signal  contained  in  an  entire  seismogram  and 
that  we  completely  bypass  the  problem  associated  with  overlapping  energy  of  wavegroups  be¬ 
longing  to  different  orbits.  It  is  not  necessary  that  the  analyzed  seismogram  should  begin  with 
the  first  minor  arc  arrival  Rj.  It  may  begin  at  any  time  tQ,  providing  that  the  interval  t}  to  tQ 
is  sufficiently  long;  190  min.,  the  circumferential  travel  time  of  the  slowest  mantle  Rayleigh 
waves,  is  a  conservatively  safe  estimate. 

The  entire  procedure  has  been  tested  first  on  synthetic  data;  application  to  the  analysis  of 
observed  seismograms  is  described  in  Sec.  C  below.  The  synthetic  seismogram  shown  in 
Fig.  Ill  —  3  has  been  obtained  by  superposition  of  the  fundamental  spheroidal  modes  computed  for 
an  Earth  model  perturbed  by  the  amounts  listed  in  Table  III— 1  with  respert  to  the  reference 
model.  The  objective  of  this  test  is  to  find  out  how  well  our  procedure  will  perform  in  deter¬ 
mining  the  values  of  these  perturbations  if  we  use  the  transfer  function  T  for  the  reference 
model  in  the  starting  iteration.  The  test  is  a  rather  difficult  one,  as  the  perturbations  in  the 
shear  velocity  offset  each  other. 

The  top  trace  in  Fig.  Ill -4(a)  is  the  "observed"  seismogram  displayed  for  the  time  interval 
corresponding  to  the  box  in  Fig.  III-3.  The  next  trace  is  the  result  of  convolution  of  the  com¬ 
plete  seismogram  of  Fig.  Ill  —  3  with  the  transfer  function  for  the  starting  model.  Even  visual 
inspection  reveals  that  there  are  substantial  differences  between  these  two  time  functions.  Their 
values  are  displayed  in  the  third  trace  and,  indeed,  the  amplitudes  in  the  difference  trace  are 
comparable  to  those  of  the  signal.  Despite  the  magnitude  of  differences,  the  procedure  con¬ 
verges,  and  the  bottom  trace  shows  that  the  residual  function 

[S(t)  -  S(t)  *  T(t)  I, 

t,  2 
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TABU  lll-l 

PERTURBATION1*  TO  THE  UURfNCl  A'ODtL  UUD 
IN  THF  SYNTHETIC  TEST T 


Depth  Range 
ft-  ff ) 

1  9  to  80 

i 

-  +■  - 

s 

11  ft  if 

0  ’ 

t  • 

'  q 

u 

-0  (J030 

SO  to  220 

j 

0  1 

i 

0  0025 

220  to  400 

0.  1 

1 

-U  0032 

400  to  670 

, 

-0  1 

i 

-0  0022 

_ L 


i  T ne  Eai  lii  model  uied  to  geneiate  the  .yntheric  seiirr,ogtarr 
sh own  in  E  ig  111-3  was  obtained  by  peiturbing  the  sbeat  ve¬ 
locities  and  attenuation  of  model  PREM  of  CteiewonsTi  and 
Andetson^  by  the  amounts  listed  obove.  Through  the  I  in  — 
eat  ized  inversion  of  the  seismogram  in  Fig.  111-3  we  shall 
attempt  to  recover  these  values,  using  the  transfer  function 
computed  for  PREM  in  the  starting  iteration.  The  quantity 
q  is  the  reciprocal  of  the  quality  factor  Q, 


is  very  .small  and,  therefore,  our  transfer  function  obtained  as  the  result  of  inversion  should 
represent  the  "true"  dispersion  and  attenuation  rather  w  ell.  The  frequency-domain  equivalent 
of  Kit;.  Ul-T(a)  is  shown  in  Kip.  Ill— lib),  where  t he  moduli  of  the  spectral  amplitudes  aredisplayed. 

Throughout  all  iterations  that  led  to  the  final  result,  the  differential  kernels  (Mr,  w)  were 
kept  unchanged,  even  though  they  arc  functionals  of  the  model  and  the  model  ehanged  during  in¬ 
version.  It  is  desirable  to  recompute  the  eigenfrequeneies  and  differential  kernels  after  the 
first  pass,  and  repeat  the  entire  process,  figure  I  i  I  -  A  demonstrates  that  the  effect  of  the  non¬ 
linearity  is  discernible.  The  starting  transter  function  is  now  computed  for  the  final  model  of 
"1'ass  1."  Note  that  amplitudes  of  the  difference  at  the  beginning  of  Pass  Z  are  greater  than  at 
the  end  of  Pass  t.  This  means  that  the  effect  of  nonlinearity  is  important  if  one  considers  deri¬ 
vation  of  the  Karth  model  to  he  the  principal  objective  of  the  procedure.  If  only  measurement 
of  dispersion  and  attenuation  is  of  importance,  then  the  accuracy  of  the  result  obtained  at  the 
end  of  Pass  1  will  be  sufficient  for  practical  purposes. 

Figure  ITI-Mal  shows  comparison  of  the  phase  velocities  measured  after  Passes  1  and  3 
with  the  exact  results;  the  largest  'inference  after  Pass  1  is  about  -  parts  m  100,000,  and  after 
Pass  2  the  results  agree  t<>  five  significant  figures.  The  differences  in  attenuation,  shown  in 
b  ig.  III-blbl,  are  more  substantial  and  for  Pas.-,  i  they  reach  Z  percent  at  b00  s.  One  can  ex¬ 
pect,  however,  that  the  magnitude  of  this  error  caused  by  the  nonlinearity  of  the  procedure  is 
small  iri  comparison  with  other  sources  of  error  that  would  be  encountered  jn  processing  actual 
data. 

In  Kig.  Ill-7(a-b)  we  compare  the  retrieved  structures  with  the  known  perturbation.  The 
velocities  after  Pass  1  are  very  close  to  the  correct  answer  -  the  largest  deviation  is  0.01  km's 
or  10  percent  of  the  perturbation.  However,  there  are  major  differences  in  the  (.J  structure, 
m  particular,  the  sign  of  the  perturbation  in  the  lid  is  opposite  to  the  correct  one.  The  structure 
derived  m  Pass  1  is  extremely  close  to  the  true  structure. 
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with  tin  e;  t tie  gradual  decrease  in  tile  frequency  content  suggests  that  it  represents  an 
attenuated  seismic  signal,  must  likely  associated  with  the  energy  of  overtones.  The  equivalent 
trace  in  Kip.  lll-'Ub),  representing  the  power  spectrum  of  the  residual  trace,  indicates  that  no 
correlation  with  the  spectral  peaks  of  the  fundamental  spheroidal  mode  remains.  We  have  dem¬ 
onstrated,  therefore,  that  the  method  does  indeed  succeed  in  properly  modifying  the  transfer 
function  through  perturbation  of  the  structural  parameters. 

The  purpose  of  1  ig.  111-10  is  to  demonstrate  that  dispersion  measured  in  such  a  way 
is  in  general  agreement  with  the  results  obtained  using  standard  techniques.  The  con¬ 
tinuous  lines  mi  that  figure  represent  results  of  phase-velocity  measurements  obtained  from 
the  phase-delay  function  of  the  appropriately  selected  fragment  of  the  autocorrelogram  of  the 

1  i 

observed  seismogram  the  smooth  broken  lines  are  the  values  resulting  from  application  of 
the  technique  scribed  in  this  report.  It  is  clear  that  there  is  an  overall  correspondence  be¬ 
tween  thos>  .o  sets  of  answers;  the  smooth  linos  represent  the  trend  of  the  results  from 
the  phase  d-iav  analysis,  but,  in  addition  to  being  physically  realistic,  they  also  extend  to  sub¬ 
stantially  longer  periods. 

figure  111-1  1  demonstrates  the  repeatability  of  the  results.  Two  strong  earthquakes  in  the 

Solomon  Islands  IM  6.1  and  7.1)  occurred  within  22  h  of  each  other;  their  locations  were  ex- 
s 

tremoly  close,  and  visual  inspection  of  several  recordings  indicates  that  their  source  mech¬ 
anisms  were  also  similar.  The  dispersion  curves  shown  in  Fig.  ill -9 (a)  are  identical  for  a 
large  range  of  angular  order  numbers,  and  the  maximum  difference  at  short  periods  is  of  the 
order  of  0.0  1  percent.  This  agreement  is  even  more  remarkable,  as  the  delay  with  respect  to 
the  origin  time  was  different  for  the  two  earthquakes  because  of  the  difference  in  their  magnitude. 

Bv  a  remarkable  coincidence,  we  have  found  recordings  of  three  earthquakes  on  nearly  the 
same  great-circle  path  with  a  station  at  Nana,  Bern  (see  Fig.  111-12);  the  poles  of  the  individual 
paths  are  within  100  km  of  each  other,  Unlike  in  the  case  shown  in  Fig.  IH-11,  the  locations  of 
the  events  were  quite  different.  We  do  not  know  the  source  mechanisms  of  these  earthquakes, 
but  it  is  highly  unlikely  that  they  would  be  similar.  Although  there  is  an  overall  similarity  be¬ 
tween  the  three  dispersion  curves,  the  differences  are  several  times  greater  than  the  maximum 
difference  in  Fig.  HI- 11.  This  effect  can  be  explained  by  the  results  of  Woodhouse,  presented  in 
Sec.  II  below.  It  shows  that  the  perturbation  in  the  period  of  a  particular  normal  mode  depends, 
in  general,  on  the  location  of  the  source  and  receiver  as  well  as  on  the  source  mechanism. 

We  have  obtained  results  for  37  great-circle  paths  and  performed  decomposition  of  the  ob¬ 
served  phase  velocities  and  attenuation  in  terms  of  the  geometrical  "pure  path"  analysis  first 
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proposed  by  Toksdz  and  Anderson,  and  later  used  by  Kanamori  and  Dziewonski  as  well  as 

several  other  authors. 

•Stable  results  were  obtained  to  periods  as  long  as  600  s,  roughly  twice  the  range  achieved 
in  the  earlier  analyses.  However,  it  appears  from  the  work  of  Woodhouse  that  there  is  a  dis¬ 
tinct  possibility  that  the  geometrical  ray  approximation  may  lead  to  serious  errors,  particu¬ 
larly  at  very  long  periods  (wavelengths).  For  this  reason,  we  do  not  think  that  publication  of 
our  current  results  would  be  particularly  useful. 

A.  M.  Dziewonski 
J.  M.  Steimt 
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!>.  A  NOV-ASYMPTOTIC  "KAY  THEORY"  FOK  SlRFACi;  WAVKS 

Ml)  FREE  OSt'll. RATIONS 

l  ln-  cifi'i't  of  lateral  heterogeneities  on  the  dispersion  of  surface  waves  is  known  to  be  sub- 

12  -  i  5 

stantial  in  tin-  Earth.  and  has  been  the  subject  of  many  studies.  Tile  interpretation  of 

such  dispersion  data  invariably  lias  been  performed  under  the  assumption  that  the  dispersion 
characteristics  are  those  appropriate  for  a  model  with  the  structure  which  is  the  average  of 
till'  true  structure  along  the  great-circle  path  for  which  the  data  were  obtained.  For  a  given 
regionalization  of  the  llarth,  this  ray  theoretic  approximation  enables  one  to  invert  for  the  dis¬ 
persion  characteristics  of  each  region,  assumed  uniform,  and  thus  to  quantify  the  variations 
in  structure  among  the  regions.  We  shall  refer  to  this  approach  as  the  "pure  path"  technique. 

Another  approach  to  the  treatment  of  the  effects  of  lateral  heterogeneity  is  through  the 
theory  of  spectral  splitting  of  free  oscillations.1^  We  shall  refer  to  this  as  the  "splitting 

technique,"  though  it  is  a  technique  which  has  not  been  applied  to  actual  data. 

1  8 

It  has  been  shown  uy  .Iordan  that  a  convenient  connection  between  the  two  techniques  may 
be  established  through  the  spectral  "location  parameter"  A.,  namely  the  difference  in  frequency 

between  the  centroid  of  a  spectral  peak  corresponding  to  a  given  (fundamental)  normal  mode  and 

1  8 

the  corresponding  eigenfrequency  of  a  reference  Earth  model,  in  particular,  it  has  been  shown 
that  asymptotically,  in  the  limit  of  large  angular  order,  the  location  parameter  for  a  particular 
source  and  receiver  is  simply  the  average  along  the  appropriate  great-circle  path  of  6  a  jocap 
which  is  defined  at  each  point  of  the  globe  (O,  A)  as  the  frequency  shift  appropriate  for  a  spheri¬ 
cal  model  possessing  the  structure  locally  beneath  that  point  (the  angles  O,  tf>  will  be  used  to  de¬ 
note  colatitude  and  longitude,  respectively).  This  result  provides  the  theoretical  justification 
for  the  pure  path  technique,  and  is  intuitively  very  appealing.  In  order  to  obtain  it,  however, 
severe  approximations  must  be  made.  In  particular,  it  must  be  assumed  that  the  wavelengths 
characteristic  of  the  lateral  variations  in  structure  are  much  longer  than  those  of  the  surface- 
wave  mode  under  consideration,  a  situation  which  probably  does  not  obtain  in  the  Earth;  indeed, 
a  remarkable  and  graphic  demonstration  of  the  inadequacy  of  the  ray  theoretic  approximation, 
using  data  from  the'  IDA  network,  is  presented  by  Dziewonski  and  Steim  in  Sec.  C  above. 

It  is  the  purpose  of  this  report  to  show  that  the  result  from  the  theory  of  splitting  for  the 
spectral  location  parameter  may  be  recast,  without  any  further  assumptions,  into  a  form  sim¬ 
ilar  to  that  used  in  the  pure  path  technique.  The  difference  is  that,  instead  of  depending  only 
upon  6a  \  depends  upon  three  local  functionals  of  Fiarth  structure  |6a(’\o,  <p),  i  =  1,3) 

and,  in  addition,  the  path  average  must  be  performed  using  not  the  geometrical  path  lengths  in 
each  region  but  three  other  parameters  which  are  fairly  readily  calculated,  and  for  which  for- 
r-iln.i  will  be  given.  Thus,  it  is  possible  to  quantify  the  effect  on  the  location  parameter  of 
regional  structural  variations  in  a  manner  which  is  similar  to  that  assumed  in  the  pure  path 
technique,  but  which  does  not  involve  the  approximations  entailed  by  that  technique. 

Hi  shall  confine  our  attention  here  to  a  particular  free  oscillation  multiplet  nS^  or 
when-  n  is  radial  order,  /  is  angular  order,  and  S  and  T  refer  to  spheroidal  and  toroidal 

■  "dr.-,,  equivalently,  we  shall  consider  a  Rayleigh  or  I.ove  wave  at  some  particular  wavelength. 

1  8 

From  the  theory  of  splitting,  it  may  be  shown  that  for  a  given  source-receiver  pair  the 
location  parameter  is  given  by 
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whore  the  summations  are  over  azimuthal  order  for  -I  <  m,  in'<  l  .  In  Eq.  (111-11),  Rm  will 
he  referred  to  as  the  receiver  vector  and  is  given  by 
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and  S  will  be  referred  to  as  the  source  vector,  and  is  given  by 
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[n  Eqs,  (III- 12)  and  (111-13),  \  (0,0)  are  the  generalized  spherical  harmonics  of  Phinney  and 

1  q  r 

Burridge,  and  O  ,  <i  ,  O  .  and  o  are  the  colatitude  and  longitude  of  receiver  and  source. 

a  R^fO,  0),  S„  =  SN(0,  0)  depend  upon  the  orientation  of  the  component  measured  by  the  seis¬ 
mometer  and  the  moment  tensor  of  the  source,  respectively,  and  are  given  in  Table  III- i.  In 
Eq.  (Ill- 11),  Hm,m  are  components  of  the  splitting  matrix.16,1  7  The  splitting  matrix  will  in¬ 
volve  components  dependent  upon  the  Earth's  rotation  and  ellipticity  H^r,^,  together  with  those 

relating  to  the  difference  between  the  reference  Earth  model  and  the  actual  Earth.  H  is 

m' m 

readily  calculated,  and  here  we  shall  be  concerned  with  only  the  remaining  terms.  Therefore, 
we  shall  write 


II  ,  =  Il,r,e)  +  H  , 

m  m  m'm  m  m 


Let  us  suppose  now  that  the  Earth  differs  from  the  reference  model  by  amounts  6p,  ok  ,  6g, 
O0Q,  ‘Sdpd’g  in  density,  bulk  modulus,  shear  modulus,  gravitational  potential,  and  its  radial  de¬ 
rivative,  respectively,  and  let  us  denote  these  model  perturbations  by  the  vector  6m  =  6m(r,  9,  0) 
The  perturbations  may  fce  expanded  in  spherical  harmonics: 


6m  =  \  fim1:Y1(9,0) 

—  — '  — s  s 

st 

17  —  t 

Woodhouse  and  Dahlen  have  given  expressions  for  in  terms  of  6m  (r),  which  we  need 

not  repeat  here  except  to  note  that  they  take  the  following  form: 


»mm.  =  Yj  ysmrr,,t  \  dm^(r)  ’  {M<0)< >')  t  s(s  i  1)  M(1)(r) 
st 


+  [s(s  4  1112  M<2)<>-)l  f2dr 


(111-14) 
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TABLE  HI-2 

COEFFICIENTS  FOR  EVALUATING  THE  SPECTRAL  LOCATION  PARAMETER 


N 

SN  a  SN(0,0) 

0 

k0^r 

VarUMrr  +  '/2F(M0e  +  M^)J 

tl 

h  |  (V  ±  iW)  (TVg  -  iv^) 

k  (X  ±  IZ)  (TM  _  —  iM  ) 

1  rd  r<J> 

±2 

" 

k2r;'(V  ±  iW)|(Meg-  M^)  ±  2iMflf) 

U,  V,  W  are  scalar  eigenfunctions^ 
evaluated  at  the  Earth’s  surface 

U,  V,  W  are  scalar  eigenfunctions^ 
evaluated  at  the  source  radius  r 

s 

,  =  J.,2 t  +  I  (I  *  n).',l/2 

n  2n  4ir  (I  —  n)l 


F  =  r"'(2U  -  F(f  +  1)  V! 
X  =  V  +  r"'(U  -  V) 

Z  =  W  -  r-1W 


NOTES:  M  ,  M^g,  etc.  are  components  of  the  moment  tensor. 

v  ,  v_,  v  are  components  of  the  "instrument  vector." 
r  0  <p 

Viz.  a  unit  vector  in  the  direction  of  motion  sensed  by 
the  instrument,  multiplied  by  the  instrument  response. 

Dot  denotes  differentiation  with  respect  to  radius  r. 


where  M^'(r)  are  given  by  known  expressions  in  terms  of  scalar  eigenfunctions  and,  what  is 
important  here,  are  independent  of  s.  The  coefficients  y™™  1  are 

mm't  _  f  m*  tym' 

s  J  t  st 

where  the  integration  is  over  the  unit  sphere. 

Now,  ifY]  denotes  the  gradient  operator  on  the  sphere: 

A  A 

—  1  =®9q  +  cosec9  _ 

we  may  use  the  divergence  theorem  on  the  .  sphere  to  write 

Z  [ s ( s  +  l))1  T^^'^m^r)  =  Z  6m!.lr)  §  Y^-V2)*  (Y™*  Yfm')  d« 

st  st 

=  C  am(r,  0,  0)  k.nim'(e,  0)  dft 


with 


kmrn  (O,  0)  =  (-V2)1 


and  Kq.  (ill-14)  becomes 
2 

H  ,  =  \  {  6m(r,  O,  0)  k.mm  (O,  0)  r2drdSl 

i=0 

and  the  location  parameter  is  given  by 
2 

A  =  X  ^  6w(l)(0,  0)  KjfO,  0)  d!2 
i=0 


with 


K.(0,  0) 


2  Rm'sm*km’m<0,0) 

mm  ' _ * _ 

[2i(<  +1)1'  2  RmSm* 


m 


(111-15) 


6cc(l)(0,0)  5{2 HI  +  1))M  6m(r,0,0)  •  M(l)(r)  r2dr  . 

The  factors  [2f((  +  1 )  1 1  have  been  introduced  to  normalize  the  kernels  in  a  convenient,  but  some¬ 
what  arbitrary,  way. 

Equation  (111-15)  is  the  principal  theoretical  result  of  this  report,  and  states  that  A  is  given 
by  the  weighted  average  over  the  surface  of  the  Earth  of  three  local  functionals  of  Earth  struc¬ 
ture,  6cu^,  with  respect  to  kernels  which  are  readily  calculated,  provided  source  and  receiver 
are  known.  is  identical  to  $«iocai  of  Jordan.18 

The  kernels  1^(0,  0)  are  most  easily  evaluated  by  making  use  of  a  coordinate  system  in 
which  the  point  (0,  0)  is  located  at  the  pole.  It  may  then  be  shown  that 


Kt(e,0)  =  £  (  — 


+  1, 


4ir 


M<.  =>  p 
ai 


NPf  <V 


,  NM, 


A  ) 
Pr 


NMN1 

x  exp  {i  [N0pr  -  N'0ps  +  M(0sp  -  0 rp) ) } 

M  NM  ‘ 

where  a.  are  given  in  Table  III  — 3,  and  are  generalized  associated  Legendre  functions. 

The  angles  A,  0  are: 

Aps  =  angular  distance  of  (0,  0)  from  (0g,  0g) 

A  =  angular  distance  of  (0.  0)  from  (0,0) 
pr  0  r'  r 

0pg  =  azimuth  of  (0,  0)  from  (0g,  0g) 

0  =  azimuth  of  (0,  0)  from  (0,0) 

pr  r  r 

0gp  =  azimuth  of  (0g,  0g)  from  (0, 0) 

0rp  =  azimuth  of  (0  0^)  from  (0 , 0). 
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TABLE  111-3 

COEFFICIENTS  FOR  EVALUATING  THE  PATH  FRACTION  KERNELS 


M 

a. 

i 

M  =  0 

M  =  ±1 

M  =  ±2 

;  =  o 

1 

- 

- 

i  =  i 

1 

-1/2 

- 

;  =  2 

3/2 

I2t(t  +  l)]'1  -  1 

1/4  -  (2f (/  +  |)f' 

Since  the  summations  over  M,  N1  are  only  from  -2  to  >2,  and  those  over  N  from  -1  to  -1,  the 
kernels  can  be  rapidly  computed.  Indeed,  for  a  vertical  instrument  and  an  explosive  source, 

KQ  is  simply  proportional  P^lcosA^)  P^(cos  where  P{  are  Legendre  polynomials. 

Some  examples  of  these  kernels  KQ (Q ,  0)  are  depicted  in  Figs.  Ill  —  1 3 ( a)  through  (c).  These 

figures  represent  the  surface  z  =  Kq(0,  0)  plotted  over  a  Mercator  projection  of  the  Earth's 

surface,  in  which  a  source  and  a  vertical  instrument  lie  on  the  equator,  a  distance  of  108°  apart. 

The  source  is  located  at  the  rightmost  (i.e.,  eastmost)  end  of  the  plots;  the  location  of  the  source 

and  receiver  and  their  antipodes  are  evidenced  by  the  peaks  in  Figs.  III-l 3(a)  and  (c).  For  QS45 

[Fig.  III-l  3(c)  shows  the  kernel  for  an  explosive  source],  it  is  clear  that  the  ray  theoretic  re- 
1 8 

suit  of  Jordan  is  emerging;  the  surface  wave  averages  local  structure  along  the  great-circle 
path  (in  this  case,  the  equator),  since  over  the  rest  of  the  Earth  the  rapid  oscillations  will  tend 
to  cancel  in  the  integral.  It  is  notable,  however,  that  the  kernel  is  still  very  broad,  and  does 
not  become  uniformly  small  toward  the  poles,  as  might  have  been  expected  on  the  basis  of  the 
asymptotic  result.  At  lower  frequencies  it  is  clear  that  very  significant  deviations  from  ray 
theory  can  occur,  and  that  they  will  depend  upon  source  and  receiver  locations,  and  also  upon 
the  orientation  of  the  instrument  and  upon  the  source  mechanism.  This  last  is  illustrated  by 
Figs.  III-l  3(a)  and  (b),  in  which  the  only  difference  is  that  in  the  first  the  source  is  explosive, 
whereas  in  the  second  it  is  vertical  strike  slip  faulting  at  an  angle  45°  to  the  equator. 

Let  us  now  suppose  that  the  heterogeneous  Earth  comprises  J  laterally  homogeneous  re¬ 
gions  FL.  Then  we  have 

*  =  £  £  Pj<l)<5“j(l)  (HI-16) 

i=0  j=l 


where  6u>^  are  those  appropriate  to  the  j1^1  region  and 


K.le- 


0)  dll 


R. 

1 


We  may  refer  to  p,^  as  the  path  fraction  of  the  i^  kind  for  the  j1*1  region.  Indeed,  if  the  path 

J  4L 

fraction  of  the  zeroth  kind  is  replaced  by  the  geometrical  path  fraction  for  the  j  n  region,  and 
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the  path  fractions  of  the  first  and  second  kinds  are  neglected,  Eq.  (111-16)  reduces  to  Jordan's 
asymptotic  result  -  the  pure  path  approximation.  The  significance  of  10 q.  (111-16)  is  that  the 
data  may  be  interpreted  in  a  way  exactly  analogous  to  the  pure  path  technique,  without  making 
any  approximations  beyond  those  inherent  in  perturbation  theory. 

Let  us  define  the  functions  m(0,  <S)  to  be  unity  in  the  j1*1  region,  and  zero  elsewhere.  If 
these  functions  are  expanded  in  spherical  harmonics: 


18 


n.(0,  0)  =  E  n  *  Y  t<0. 
J  ^  ,1^  s 

st 


We  obtain  the  following  expression  for  the  path  fractions: 

.  (i)  1  V  ,  cm*/^  ,  f  s(s  +  1)  |i 

Pj  *  p  l  R  0  r*  ^  r  S  <6s-  *s>  2i  (f  +  l)1 


X(-l)m  (2  t 


,  , .  ,  2s_t  1  .  l/2 /f  l  s\  /  t  I  s  \  n 
ir  (o  0  0/(-m'  m  m-m')  js 


where  Wigner  3-j  symbols21  have  been  introduced.  While  the  evaluation  of  the  quai  tities  p/11 
is  somewhat  intricate,  it  is  not  expensive  computationally  and  we  have  developed  algorithms 
for  their  calculation.  Work  proceeds  on  using  these  parameters  in  the  interpretation  of  dis¬ 
persion  data.  j.h.  Woodhouse 

T.  P.  Girnius'f 
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Fig.  Ill  —  1 .  Rayleigh- and  Love-wave  radiation 
patterns  for  various  combinations  of  explosive 
and  thrust  faulting,  expressed  as  RfB.  the 
ratio  of  seismic  moment  of  faulting  component 
to  that  of  explosion.  Shaded  portions  denote 
source  phase  of  -II,  unshaded  a  phase  of  zero. 
Dashed  circle  indicates  amplitudes  from  ex¬ 
plosion  only.  Shape  of  patterns  remains  un¬ 
changed,  increasing  in  amplitude,  as  Rp-jj 
exceeds  3.0. 
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Fig.  III-2.  Minimum,  maximum,  and  azimuthally  averaged  Rayleigh-wave 
spectral  amplitudes  as  a  function  of  Rp.^. 


0  2  4  6  a  10 

TIME  FROM  START  OF  RECORDING  (h) 


Fig.  1 1 1  -  3 .  Synthetic  seismogram  for  "unknown"  structure.  Seismogram  contains 
modes  of  from  oS2  to  0S97;  source  was  at  depth  of  70  km  and  receiver  at  distance 
of  60°  from  epicenter.  Theoretical  seismogram  has  been  convolved  with  LP 
WWSSN  response  and  a  tapered  low -pass  filter  applied  to  remove  energy  at  periods 
shorter  than  100  s.  Using  technique  described  in  this  report,  we  shall  derive  ele¬ 
ments  of  elastic  and  anelastic  Earth  structure  that  are  consistent  with  information 
on  dispersion  and  attenuation  of  surface  waves  contained  in  this  recording. 
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Fig.  Ill- 4(a).  Top  trace  is  segment  of  syn¬ 
thetic  seismogram  shown  in  Fig.  Ill- 3  within 
area  outlined  by  box.  Next  trace  is  trun¬ 
cated  result  of  convolution  of  seismogram 
in  Fig.  Ill-  3  with  transfer  function  T  for 
starting  Earth  model.  If  transfer  function 
were  consistent  with  "true"  dispersion  and 
attenuation,  then  third  trace  would  be  zero; 
it  clearly  is  not.  However,  we  improve 
transfer  function  in  iterative  procedure, 
such  that  amplitude  of  difference  trace 
shown  at  bottom  of  figure  is  practically 
zero. 
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Fig.  Ill- 4(b).  Frequency-domain  rep¬ 
resentation  of  information  shown  in 
Fig.  III-4(a).  Only  power  spectrum  is 
shown. 
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starting  iteration 
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Fig.  Ill- 5.  Same  as  Fig.  Ill  4(a)  except 
for  improved  starting  model  obtained 
during  first  stage  (pass)  of  inversion 
procedure.  Discrepancy  between  last 
trace  in  Fig.  III-4(a)  and  third  trace  in 
this  figure  is  result  of  nonlinear  re¬ 
lationship  between  functionals  of  a 
model  and  perturbations  in  structural 
parameters. 


time  after  start  of  record  (h) 


Fig.  Ill- 6.  Comparison  of  (a)  phase 
velocities  and  (b)  attenuation,  ob¬ 
tained  through  linearized  inversion 
of  waveform  data  with  exact  results. 
Values  for  starting  model  represent 
reference  level. 
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Fig.  III-7.  Comparison  with  exact 
result  of  (a)  shear  velocity  and 
(b)  shear  attenuation  profiles  re¬ 
trieved  through  linearized  inver¬ 
sion  of  waveform  data. 
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Fig.  ITT-  S .  Recording  of  Solomon  Islands  earthquake  of  29  July  1977  by  an  IDA 
instrument  at  Halifax,  Nova  Scotia.  Difference  S(t)  —  S(t)  #  T(t)  is  minimized 
within  area  outlined  by  box.  Symbols  Rn  denote  individual  arrivals  of  wave- 
groups  of  mantle  Rayleigh  waves. 
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Fig.  III-9(a).  Top  trace  is  segment  of  synthetic  seismogram  shown  in  Fig.  Ill- 8 
within  area  outlined  by  box.  Next  trace  is  truncated  result  of  convolution  of 
seismogram  in  Fig.  Ill- 8  with  transfer  function  T  for  reference  Earth  model. 
Third  trace  reflects  degree  of  misfit  of  reference  Earth  model  and  structure 
for  particular  path;  deviations  of  this  magnitude  would  be  observed  in  source 
mechanism  studies,  if  excitation  kernels  are  path  independent.  After  struc¬ 
ture  is  "corrected,"  amplitude  of  difference  trace  at  bottom  of  figure  is  sig¬ 
nificantly  diminished  and  does  not  show  correlation  with  arrivals  of  individual 
wavegroups. 
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PERIOD  III 

Fig.  [[[-9(b).  Frequency-domain  representation  of  Fig.  Itl-9(a);  only  power 
spectra  are  shown.  Notice  that  spectrum  of  difference  trace  after  inversion 
appears  to  be  white  and  that  there  is  no  correlation  with  spectral  peaks  of 
fundamental  spheroidal  mode. 


Fig.  [[[-10.  Comparison  of  phase -velocity  measurements  obtained 
by  linearized  inversion  technique  (FIT)  described  in  report  (bro¬ 
ken  lines)  with  those  derived  by  traditional  phase-delay  method 
(solid  lines).  Curves  obtained  by  application  of  FIT  are  naturally 
smooth  and  can  easily  be  differentiated  to  obtain  group  velocities, 
for  example. 
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Fig.  Ill- 11.  Repeatability  of  measurements  is  demonstrated  by  results 
of  measurements  performed  on  recordings  of  two  earthquakes  of  nearly 
identical  location  and  similar  source  mechanism. 
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Fig.  Ill-  1 2.  Measurements  for  three  earthquakes  recorded  at  Nana,  Peru. 
Great-circle  paths  were  nearly  identical  (their  poles  are  within  100  km  of 
each  other),  but  epicenters  were  in  different  seismic  zones  (Tonga,  New 
Guinea,  and  Japan)  and  source  mechanisms  were  most  likely  different. 
Woodhouse  (Sec.  D,  this  issue)  shows  that  for  a  given  great-circle  path 
and  a  fixed  receiver,  observed  period  of  free  oscillation  depends,  in  gen¬ 
eral,  on  position  of  source  and  its  mechanism. 
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